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This  report  was  prepared  under  a  project  funded  by  the  Montana  Department  of  Natural 
Resources  and  Conservation.    The  Department  makes  no  warranty,  express  or  implied, 
nor  assumes  any  legal  liability  or  responsibility  for  the  accuracy,  completeness,  or 
usefulness  of  any  information  apparatus,  product,  or  process  disclosed,  nor 
represents  that  its  use  would  not  infringe  on  privately  owned  rights.  Reference 
herein  to  any  specific  commercial  product,  process,  or  service  by  trade  name, 
trademark,  manufacturer,  or  otherwise,  does  not  necessarily  constitute  or  imply  its 
endorsement,  recommendation,  or  favoring  by  the  Department  of  Natural  Resources  and 
Conservation.    The  reviews  and  opinions  of  authors  expressed  herein  do  not 
necessarily  state  or  reflect  those  of  the  Department. 


Introduction 

Four  %  of  the  total  United  States  energy  consumption  and  14  %  of  energy  use  in 
the  residential  sector  is  used  for  domestic  hot  water  production  (Hirst,  1980) . 
With  higher  prices  for  conventional  fuels  used  for  this  purpose,  the  economics  of 
the  substitution  of  solar  becomes  more  attractive.    A  computer  simulation  of  the 
cost  of  installing  solar  hot  water  (SHW)  heating  systems  would  be  useful  to  home 
owners  attempting  to  minimize  energy  related  expenses. 

previous  computer  evaluations  of  solar  domestic  hot  water  heating  such  as  TRNSYS 
(Klein  et  al.,  1976),  LASL  (Balcomb,  Hedstrom,  1976)  and  SQLCOST  (Connolly  et  al., 
1976)  have  focused  on  the  relation  of  the  collection  area  and  local  climate  to  the 
variations  in  the  percentage  (or  "solar  fraction")  of  hot  water  provided  by  solar 
systems.    Each  of  these  hourly  simulations  must  be  run  on  mainframe  computers.  This 
is  clearly  beyond  the  capabilities  of  lay  persons  who  need  such  design  tools. 

DNRC  has  developed  an  easy  to  use  solar  domestic  hot  water  program  designed  to 
run  on  microcomputers.    The  DNRC  program  is  based  on  the  FCHART  method  (Duff  ie  and 
Beckman,  1980) ,  a  widely  accepted  approach  to  estimating  the  solar  fraction  by 
simulating  the  performance  of  a  solar  system.    This  method  allows  prediction  of 
comparative  solar  vs.  conventional  fuel  costs  based  on  predicted  performance,  and 
use  of  monthly  data  inputs.    The  simulation  model  developed  by  DNRC  is  an 
interactive  BASIC  computer  program  consisting  of  four  discrete  sections:  an 
analysis  of  energy  demand,  a  collector  performance  algorithm,  a  specification  of 
economic  parameters,  and  a  series  of  sensitivity  tests.    The  availability  of  recent 
solar  data  for  30  sites  around  Montana  permits  analysis  by  geographic  region. 
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Comparison  to  Other  Programs 


A  simplified  FCHART  has  been  publicly  available  in  the  Montana  solar  ^..i 
for  over  two  years.    It  allows  analysis  of  both  active  solar  hot  water  and  space 
heating  systems  using  a  hand  calculator  and  work  sheets.    The  method  is  simplified 
for  ease  of  use,  but  may  be  cumbersome  for  some  lay  persons. 

The  DNRC  SHW  analysis  model  is  a  public  domain  program,  formatted  to  run  on 
Apple  II  microcomputers,  which  are  widely  available  at  schools,  universities  and 
businesses  around  the  state.    It  uses  the  FCHART  algorithm  (see  Appendix  A)  in  a 
detailed  form  to  predict  solar  hot  water  system  output  and  costs.    The  DNRC  model 
makes  several  improvements  to  the  simplified  FCHART: 

1)  The  program  is  interactive  in  a  friendly  fashion.    Questions  are 
displayed  on  the  screen  and  answers  are  requested  from  the  keyboard.  No 
mathematic  skills  are  required  for  the  inputs.    Computations  are  fast  and 
automatic  with  a  typical  runtime  of  less  than  five  minutes  for  a  complete 
analysis.    There  is  extensive  error  checking  in  the  program,  so  that 
erroneous  data  is  not  entered.    Default  values  are  provided  for  all 
questions  based  on  typical  responses.    It  is  possible  to  print  a  summary 
of  the  analysis  results  to  a  printer. 

2)  The  water  heating  load  is  calculated  with  good  accuracy.    Effects  of 
occupant  consumption  behavior  are  incorporated  as  well  as  a  mathematical 
model  for  predicting  effect  of  household  size  on  consumption.  Seasonal 
fluctuation  in  mains  water  temperature  has  been  studied  in  Missoula  and 
Helena,  and  an  excellent  lagged  linear  fit  to  monthly  ambient  temperatures 
has  been  derived. 
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3)  Additional  correction  factors  are  included  for  minimum  acceptable 
water  temperature  and  thermal  losses  from  the  auxilary  storage  tank,  which 
have  significant  effects  on  SHW  system  performance  (see  Appendix  A). 

4)  An  economic  analysis  section  is  included  that  allows  rapid 
optimization  of  collector  size.    Sensitivity  study  of  the  effect  of  a 
variety  of  economic  parameters  is  easily  accomplished.    These  parameters 
could  include  local  utility  rates,  energy  cost  escalation  and  general 
inflation  rates,  length  of  period  covered  by  analysis,  and  loan  interest 
rates. 

5)  Solar  and  climatological  data  is  available  for  thirty  sites  around 
Montana,  allowing  wide  application  of  the  program. 


6)  Specific  assumptions  are  made  on  the  efficiency  of  home  built  and  self 
installed  systems  based  on  measured  results  (Winn,  1981) . 


In  the  last  year  the  Montana  Power  Company  has  incorporated  FCHART  and  the 
Fowlkes  data  for  implementing  SHW  into  their  zero  interest  loan  program.    The  DNRC 
model  uses  essentially  the  same  calculation  algorithm  although  some  of  the  load 
calculations  are  more  specific.    The  program,  however,  is  designed  to  be  easy  to  use 
for  the  average  lay  person  and  is  publicly  available.    The  greatest  disparity 
between  the  two  programs  is  the  economic  assumptions  used  on  fuel  price  escalation 
rates  in  the  respective  reference  cases.    The  DNRC  assumptions  are  considerably  more 
conservative  in  this  respect,  although  the  effect  of  comparable  escalation  rates  is 
easily  studied. 
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Physical  Description  of  Simulate  solar  systems 

The  computer  simulation  example  used  in  this  report  is  based  on  the  physical 
characteristics  of  the  "drainback"  type  of  collector,  a  common  type  of  solar  hot 

i   water  system  in  Montana.    (Performance  characteristics  for  other  types  of  systems 
may  be  entered  in  the  model  if  desired.    Systems  using  air  as  the  heat  transfer 
medium  can't  be  simulated  with  DNRC's  model.)    A  schematic  of  the  system  is  shown  in 
Figure  1.    In  this  configuration,  typical  flat  plate  collectors  convert  solar 
radiation  to  heat.    That  heat  is  transferred  to  water  that  moves  via  riser  tubes  to 
a  solar  storage  tank.    This  tank  contains  a  heat  exchanger  that  heats  cold  line 
water  before  it  enters  the  conventional  hot  water  heater.    Water  can  only  reach  the 
collectors  when  the  pump  is  turned  on  by  the  thermostat.    When  the  pump  turns  off, 
the  collector  system  drains  naturally,  leaving  the  array  empty  and  f reeze-proof . 
The  water  in  the  storage  tank  is  isolated  from  tap  water,  eliminating  the  long  term 

-    calcination  problems  that  otherwise  would  occur  with  drainback  systems. 

■ 


Figure  1.  -  Drain  back  Sol  =  r  Syste.- 
Configuratio- 
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The  collectors  used  in  the  examples  here  are  liquid,  selective  surface 
double-glazed  collectors.    However,  values  for  any  liquid  collector  system  may  be 
substituted  in  the  program  by  the  user  to  allow  analysis  of  any  conventional 
configuration.    Solar  storage  and  auxiliary  tanks  are  assumed  to  be  insulated  to 
R-10.    A  ratio  of  two  gallons  of  water  for  solar  storage  is  assumed  per  square  foot 
of  collector  surface. 

Collector  Performance  Algorithm 

Thermal  performance  of  the  collectors  is  analyzed  in  order  to  determine  system 
output.    The  factors  involved  in  such  analysis  are  very  complex  thermodynamic 
processes.    The  DNRC  model  uses  the  FCHART  method  developed  by  Beckman  et  al.  (1977) 
to  simulate  collector  performance.    FCHART  is  based  on  hundreds  of  runs  of  the 
TRNSYS  simulation  and  has  been  extensively  validated  against  actual  system 
performance.    Agreement  with  experimental  results  is  considered  to  be  very  good  for 
conventional  collector  designs.    Duffie  and  Mitchell  (1982)  have  compared  annual 
performance  of  eleven  solar  hot  water  systems  with  predicted  performance  using  the 
FCHART  method.    Agreement  with  measured  solar  fractions  was  within  +  3  percent 
annually.    The  calculation  method  is  described  in  detail  in  Appendix  A. 

Data  Accuracy 

Inaccuracies  in  the  model  estimations  can  occur  at  several  points.  The 
radiation  and  climatological  data  used  in  the  calculation  procedure  has  a 
substantial  effect  on  the  results.    Inaccuracy  in  the  FCHART  collector  performance 
algorithm,  though  minimal,  may  compound  this  problem. 


The  -radiation  and  climatological  data  is  taken  from  the  Montana  Solar  Data 
Mnuai  (Fowlkes,  1981) .    The  solar  data  is  based  on  an  instrumented  thirty  site 
network  within  the  state.    Approximately  five  years  of  continuous  data  has  been 
developed  during  the  program.    Extensive  calibration  and  data  correction  efforts 
have  put  the  data  accuracy  to  within  a  +  5  percent  confidence  interval.    The  data 
recovery  rate  has  been  in  excess  of  80  percent.    The  greatest  drawback  of  the  solar 
data  is  the  relatively  short  time  period  covered. 

The  cl imatological  data  is  from  the  U.S.  Weather  Service  30-year  long  term 
averages.    A  cursory  analysis  of  this  data  for  various  Montana  locations  shows  a 
standard  deviation  of  15  percent  for  year-to-year  measured  temperatures  from  local 
long  term  averages  (NOAA,  1982) . 

Generally,  given  the  most  probable  cumulative  inaccuracies  of  the  algorithm, 
climatology  and  behavioral  aspects,  model  estimates  can  be  bounded  by  approximately 
a  ±  16  percent  confidence  interval  for  any  given  year. 

Economic  Optimization  Algorithm  -  Parameters 

The  initial  costs  of  solar  hot  water  systems  are  dependent  on  the  size  of  the 
collectors.    The  solar  fraction  increases  with  collection  area  but  at  a  steadily 
diminishing  rate.    Therefore,  it  would  be  prohibitively  expensive  to  size  a  system 
to  provide  all  household  hot  water  year-round.    Determining  the  economic  parameters 
for  the  financial  analysis  allows  estimation  of  the  most  economic  system  size. 


The  fixed  and  variable  costs  for  commercial  SHW  systems  were  derived  by 
performing  a  linear  regression  of  cost  data  compiled  by  McPherson  (1980)  and  Sawyer 
(1978) .    These  costs  were  significantly  lower  for  home  built  systems.    Some  of  the 
initial  cost  of  commercial  systems  can  be  reduced  by  purchasing  equipment  and 
performing  an  owner  installation. 

Table  1  summarizes  the  fixed  and  variable  cost  assumptions  in  the  model.  Data 
for  1978  and  1979  were  inflated  by  use  of  the  consumer  price  index  to  1982  dollars. 
Phone  calls  to  four  Montana  solar  hot  water  system  dealers  found  this  price  data  to 
be  in  close  agreement  with  current  system  costs  (Parker,  1983) . 

Table  1 

Fixed  and  Variable  Costs  of  Solar  Systems 
Cwner  Built  Commercially  Built     Commercial  Unit/ 

find  installed  and  Installed         Owner  Installed 

Fixed  $500  $1200  $700 

Variable  (ft2)  $  20  $    30  $  30 

Sources:    McPherson,  B.,  1980.  Hot  Water  From  the  Sun;  U.S  Dept.  of  Housing  and 
Urban  Development,  HUD-PDR-548,  Washington,  D.C. ; 

Sawyer,  S.W.,  1978.  Cost  Performance  and  Reliability  Patterns  of  177  Solar 
Heating  Systems;    Proceedings  of  the  Solar  Heating  and  Cooling  Operations 
Results  Conference,  SERI,  Colorado  Springs,  Colorado. 

The  cost  of  conventional  energy  was  taken  from  the  January,  1983,  residential 
rates  for  each  location  in  the  solar  data  base.    This  information  is  presented  in 
Table  2.    Although  all  gas  suppliers  have  winter  rates  for  the  first  fifteen  MCF  of 
consumption,  it  is  assumed  in  the  model  that  the  avoided  gas  cost  in  these  months 
will  be  at  the  higher  rate.    Most  customers  with  gas  water  heating  also  have  gas 
space  heat  which  will  consume  the  lower  cost  gas. 


Table  2 

Residential  Electricity  and  Natural  Gas  Prices 


Utility 


Electricity  (kWh) 


N~Gas  (MCF) 


Montana  Power  Company  $  .036 

Montana-Dakota  Utilities  .069 

Pacific  Power  and  Light1  .03/  .048  .063 

Flathead  Irrigation  Proj.2  .013/. 019 

Glacier  Electric  Coop  .028 

Fall  River  Electric  Coop  .037 

Great  Falls  Gas  NA 


$  4.672 
3.945 


NA 
NA 
NA 
NA 
5.175 


NA  -  Not  applicable 

1  Rate  Structure:    0-300  kWh  =  $.0304,  301-1000  =  $.0448,  1000+  =  $.06265,  minus 

$.00781  Exchange  Credit 

2  <900  kWh  =  $.013  kWh,  >900  =  $.019 

MDU  =  Glasgow,  Glendive,  Miles  City,  Plentywood,  Sidney 

PP&L  =  Kalispell,  Libby 

Flathead  Irrigation  Project  =  Poison 

Glacier  Electric  Cooperative  =  Browning 

Fall  River  Electric  Cooperative  =  West  Yellowstone 

MPC  =  Most  others 

Great  Falls  Gas  =  Great  Falls 

Source:  T.  Otis,  Phone  Conversation,  Public  Service  Commission,  January  4  and 
April  27,  1983. 


Twenty  years  was  assumed  as  the  most  likely  life  expectancy  for  an  active  solar 
system  of  the  type  depicted  (Beckman  et  al.,  1977).    This  agrees  closely  with  the 
"mean-time-bef ore-system  failure"  estimates  of  Goldberg  (1978) . 


Past  studies  have  underscored  the  importance  of  the  assumed  discount  rate  for 
solar-economic  models.    The  real  fuel  price  escalation  rate  is  similarly  important 
(Chandra,  1977)  .    Short-term  fuel  price  escalation  rates  were  taken  from  the 
forthcoming  Montana  Power  Company  1983  forecast  (Antonioli,  1983)  in  which  real 
rates  of  increase  were  estimated  at  6  percent  and  4  percent  annually  for  electricity 
and  natural  gas  over  a  ten-year  period.    The  long-term  real  fuel  price  escalation 
rate  is  taken  from  Nordell  (1983)  in  which  rates  for  natural  gas  and  electricity  are 


estimated  at  2.5  %  per  year.    The  Montana  Power  Company  uses  4  percent  and  6  percent 
average  real  escalation  rates  for  electricity  and  natural  gas  prices  over  the  entire 
20-year  period  in  their  analyses  of  SHW  systems  for  their  zero  interest  loan 
program.    While  such  rates  greatly  increase  solar  savings,  the  DNRC  assumptions  are 
more  conservative  for  long  term  projections.    Regardless,  these  assumptions  are  only 
little  better  than  educated  guesses. 

The  individual  discount  rate  is  input  into  the  model  based  on  whether  or  not  the 
purchaser  plans  to  borrow  money  to  pay  for  the  solar  system.    If  the  purchaser  uses 
his  own  money,  the  solar  system  can  be  treated  as  an  investment.    A  nominal  11  % 
discount  rate  is  then  assumed,  which  corresponds  to  a  best  alternative  investment 
for  an  individual  in  a  35  %  tax  bracket.    Inflation  is  assumed  to  be  8  %  per  year  in 
the  simulation.    Consequently,  the  model  uses  a  real  2.78  %  discount  rate;  that  is 
the  nominal  rate  over  inflation  (1.11/1.08  =  1.0278).    If  the  purchaser  finances  the 
system,  the  model  uses  a  13  %  interest  rate  over  the  life  of  the  investment,  or  4.6 
%  real.    This  is  lower  than  the  7.9  %  real  rate  (16.5  %  nominal)  currently  used  by 
the  federal  government  in  evaluating  energy  related  investments  (Marshall  &  Ruegg, 
1980) .    The  high  rate  used  by  the  government  tends  to  discourage  large  investments 
with  slow  rates  of  return  such  as  solar.    It  is  also  considerably  greater  than  the 
average  individual  long-term  discount  rate  (Ibbotson  and  Sinquef ield,  1979) . 

Both  state  and  federal  tax  credits  are  included  in  DNRC's  model.    The  current 
federal  tax  credit  allows  taxpayers  to  write  off  40  %  of  the  initial  investment  up 
to  $10,000.    The  state  credit  is  in  addition  to  the  federal  credit.    It  is  currently 
5  %  of  the  first  $1,000  of  the  investment  and  2  1/2  %  of  the  next  $3,000  up  to  a 
maximum  credit  of  $125  (I.R.S.,  1982;  Montana  Department  of  Revenue,  1982). 
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Operation  and  maintenance  and  pumping  and  controller  energy  costs  of  the  system 
are  set  according  to  convention  at  1  %  of  the  initial  investment  per  year  (Duff ie 
and  Beckman,  1980) .    They  are  projected  to  increase  at  the  general  rate  of 
inflation,  thus  remaining  constant  in  real  terms.    The  system  is  assumed  to  have  no 
salvage  value  after  twenty  years.    The  conventional  gas  system  and  auxiliary  gas 
backup  for  the  solar  system  are  slated  for  replacement  in  the  tenth  year. 

Economic  Optimization  -  Method 

Life-cycle  cost  analysis  is  used  to  translate  future  cost  and  benefit  streams 
into  present  value  dollars.    This  allows  direct  comparison  of  solar  versus 
conventional  options.    Optimization  of  the  solar  design  area  is  reached  at  the 
collector  area  that  maximizes  positive  differences  between  life-cycle  costs  of  the 
two  options.    This  will  result  in  the  greatest  savings  for  the  owner.    The  model 
estimates  a  variety  of  economic  measures  that  allow  assessment  of  comparative 
life-cycle  costs  and  savings,  time  to  recover  investment,  and  dollar  savings 
compared  to  dollars  invested.    All  are  presented  so  that  persons  familiar  with  a 
particular  method  will  be  able  to  compare  the  results  in  terms  of  what  they 
understand  best. 

The  present  value  of  O&M  and  replacement  costs  are  calculated  using  the  uniform 
present  worth  equation: 

P  =  A  (1  +  d) 

d  (1  +  d)n 
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P  =  A  11  +  d^"-l 
d  (1  +  d)n 

Where:  A  ■  cost  of  O&M  in  the  first  year 
d  =  discount  rate  (.0278  or  .046) 
n  =  life  of  investment  (20  years) 

The  modified  uniform  present  worth  (UPW*)  of  future  fuel  costs  uses  a  similar 
approach  modified  by  the  fuel  escalation  rate: 

10 

UPW*=Pn      z   1   + 

n=l        [(1  +  d)/(l+r)]n 

20 

I   1  

n=ll  [(l-k3)/(l+€)]n 

Where:  d  =  discount  rate  (.0278  or  .046) 
n  =  year  (1-20) 

r  =  fuel  price  escalation  rate,  years  1  to  10  (.04  gas,  .06  electricity) 

e  =  fuel  price  escalation  rate,  years  11-20  (.025) 
Using  the  nomenclature  of  Marshall  and  Ruegg,  (1980) ,  life  cycle  costs  of  the  two 
options  is  represented  by: 

LCC=P-S+M+R+E 
Where:       P  =  purchase  and  installation  costs 

S  =  salvage  value 

M  =  Present  value  of  O&M  costs 

R  =  present  value  of  auxiliary  heater  replacement 

E  =  present  value  of  energy  costs 
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Calculation  of  the  present  value  of  net  savings  results  in  a  single  number  that 
is  readily  comparable  for  ranking  purposes.    This  calculation  determines  the 
difference  between  the  life-cycle  costs  of  the  solar  investment  and  the  conventional 
alternative: 

■ 

NS  =  E  -  (P  -  S  +  M  +  R) 
Where:       E  =  reduction  in  present  value  energy  costs 

P  =  differential  purchase  and  installation  costs 

S  =  differential  salvage  value 

M  =  differential  present  O&M  value  costs 

R  =  differential  present  value  replacement  costs 

The  solar  area  that  maximizes  the  net  savings  of  the  solar  investment  is  the 
criteria  used  for  optimization.    The  program  increments  the  solar  area  by  one  square 
meter  until  the  greatest  net  savings  is  reached. 

Two  other  analytical  tools  are  used  to  rank  the  solar  system  economic 
performance.    The  savings  investment  ratio  differs  from  the  two  previous  methods  in 
that  it  provides  a  measure  of  the  economic  efficiency  (rate  of  return)  of  an 
investment.    For  positive  net  savings,  the  ratio  must  be  greater  that  one  -  the 
higher  the  ratio,  the  more  dollar  savings  realized  per  dollar  of  investment.  Using 
the  same  nomenclature  as  with  the  previous  example: 

sir  =  (e  -  m  

(P  -  S  +  R) 

The  greatest  SIR  will  often  be  that  reached  prior  to  reaching  the  most  economic 
size  of  investment,  so  the  SIR  is  not  the  preferred  method  for  optimization. 
However,  it  is  valuable  for  comparative  project  ranking. 
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The  discounted  payback  measures  the  tine  between  the  initial  point  of  investment 
and  a  time  in  the  future  when  net  accumulated  savings  offset  the  initial  costs.  The 


method  is  most  often  used  with  investments  where  a  rapid  turnover  of  investment 
funds  is  desirable.    However,  a  shorter  payback  does  not  indicate  the  most 
economically  efficient  investment. 


where: 

C±  =  initial  costs  minus  tax  credits 

Si  =  first  year  solar  savings  minus  0  &  M  costs 

d      =  discount  rate 

e      =  fuel  price  escalation  rate 

Sample  Case  parameters 

The  manner  in  which  DNRC's  program  works  can  be  demonstrated  with  a  sample 
case.    In  the  case  shown  below,  the  economics  of  solar  domestic  hot  water  were 
calculated  for  a  family  of  four  in  Helena.    The  exact  method  for  determining  the  hot 
water  load  is  detailed  in  Appendix  B.    Briefly,  the  thermal  load  varies  monthly  from 
a  low  of  1.7  x  106  Btu  in  August  to  2.1  x  106  Btu  in  March.    The  model  estimates 
that  73  gallons  of  hot  water  is  used  each  day  requiring  6,912  kwh  of  annual 
electricity  to  meet  the  demand,  including  standby  losses.    The  set  temperature  of 
the  water  tank  is  140°F.    It  is  assumed  that  external  insulation  (R-10  wrap)  has 
already  been  applied  to  the  auxiliary  tank  before  solar  is  considered. 


DPB  = 


loge  [  1  +  d  *  (e-d)] 

si  

log/l+_£\ 
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Location  of  the  residence  determines  the  amount  of  solar  energy  available  to 
meet  the  demand  and  the  cost  of  competing  conventional  fuels.    The  available  average 
daily  incident  radiation  per  square  foot  in  Helena,  Montana,  based  on  actual 
measurement  (Fowlkes,  1981),  is  shown  in  Table  3,    The  house  is  assumed  to  have  an 
unobstructed  south  facing  exposure. 

Table  3. 
Available  Insolation 
Helena,  Montana*  (Btu/ft2/day) 

Winter  Soring 
Dec.  Jan.  Feb.  Mar.  April  May 

679  805  1035  1327  1508  1362 

Summer  Fall 
June  July  Aug.  Sept.        Oct.  Nov. 

1605  1627  1640  1649  1337  964 

*   For  a  solar  collector  facing  due  solar  south,  inclined  45°  from  horizontal. 

Source:       Fowlkes,  Charless,  1982.    Montana  Solar  Data  Manual.  DNRC,  Helena, 
Montana,  p.  12  and  13. 
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Peaults  and  Sensitivity 


The  sample  computer  run  for  the  electricity  case  is  contained  along  with  the 
program  itself  in  Appendix  C.    The  optimal  collector  size  for  the  system  is  9  m2  or 
about  96  ft2  corresponding  to  the  size  necessary  to  heat  most  of  summer  hot  water 
and  70  %  of  annual  hot  water  demands.    About  2,000  kWh  of  electricity  is  used 
annually  to  provide  the  other  30  %  of  the  hot  water  demand.    The  estimated  initial 
capital  cost  of  the  system  is  $4,124,  or  $2,349  after  federal  and  state  tax 
credits.    The  life-cycle  cost  of  the  solar  system  is  $4,541  versus  $5,376  for  the 
conventional  system,  for  a  net  savings  of  $835.    The  savings  investment  ratio  is 
1.35;  the  discounted  payback  is  16  years. 


Table  4  tests  the  model  sensitivity  to  number  in  the  household  (1-6  persons) , 
hot  water  consumption  habits,  and  tank  temperature  set  point.    It  is  assumed  that 
"miserly"  users  will  consume  about  one-third  fewer  gallons  of  hot  water  than  average 
users;  "lavish"  users  will  consume  30  percent  more. 


Table  4.  -  Sensitivity  to  Number  in  Household/Habits 
Natural  gas  Electric 


Optimal 

Net 

Optimal 

Net 

ft? 

Savings 

ft2 

Savir 

#  in  Household 

4  ref.  case 

64 

.53 

$-241 

yo 

.70 

$835 

1 

21 

.44 

-752 

32 

.63 

-434 

2 

32 

.44 

-576 

53 

.66 

21 

3 

53 

.53 

-406 

75 

.68 

441 

5 

75 

.53 

-  85 

107 

.68 

1203 

6 

86 

.53 

67 

129 

.70 

1558 

4  miserly  use 

43 

.49 

-474 

64 

.66 

272 

4  lavish  use 

86 

.55 

-  5 

118 

.69 

1384 

4  120°  set  temp 

53 

.57 

-303 

75 

.73 

590 

F  =  Annual  fraction  of  hot  water  provided. 
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A  three-person  household  is  the  economic  break  even  point  for  electricity,  while 
household  habits  that  increase  use  of  hot  water  (lavish  use)  or  use  of  higher 
temperature  water  tend  to  favor  solar  investment.    A  conversion  from  natural  gas  to 
solar  is  not  cost  effective  in  most  cases.    Note  that  a  collector  of  the  same  size 
will  not  provide  the  same  solar  fraction  for  electric  and  natural  gas  systems.  This 
is  because  the  natural  gas  water  heater  cannot  be  insulated  as  well  as  the  electric 
one,  presenting  a  greater  thermal  load.    (See  Appendix  B.)    Reducing  the  water 
heater  set  temperature  to  120°  results  in  substantial  present  value  savings  ($230 
for  electricity)  at  no  investment  cost. 

Table  5  studies  variations  in  performance  and  economics  from  various  collector 
types.    The  collector  used  as  the  reference  case  is  a  good  double  glazed  collector. 
Three  other  configurations  are  modeled,  including  average  double  and  single  glazed 
collectors  and  an  evacuated  tube  model.    The  average  values  are  based  on  a  random 
sample  of  thirty  single  and  double  glazed  models  from  the  Solar  Products 
Specifications  Guide.    (Solar  Age,  1981) . 

Table  5.  -  Sensitivity  to  Collector  Type 


Natural  gas 

Electric 

Optimal 

ipi 

Net 

ft2 

ipi 

Net 

ft.2 

Savings 

Savings 

Reference  Case  64 

.53 

$-242 

96 

.70 

$  835 

Housebuilt  86 

.51 

374 

107 

.59 

1,326 

Ave.  Double  glazed  53 

.40 

-480 

86 

.58 

401 

Ave.  Single  glazed  53 

.40 

-462 

86 

.58 

406 

Evacuated  Tube  75 

.60 

-219 

107 

.78 

994 

Uninsulated  Tank  75 

.59 

123 

96 

.70 

1,283 

FrUL 

FRT  a 

Notes : 

Ref . 

Case 

.68 

.72 

Ave. 

Double  Glazed 

.81 

.67 

Ave. 

Single  Glazed 

.95 

.72 

As  seen  from  this  analysis,  there  is  little  difference  between  average  single 
and  double  glazed  collector  performance,  the  primary  deciding  factor  being  one  of 
price.    Considerably  greater  life-cycle  savings  are  possible  if  the  solar  system  is 
self-built  even  though  such  systems  are  assumed  to  have  reduced  efficiency.  The 
approximate  8  %  increase  in  efficiency  possible  from  adding  hot  water  tank 
insulation  significantly  increases  solar  savings  in  cases  where  this  will  not  be 
undertaken  without  solar  retrofit.    Evacuated  tube  collector  designs  offer  greater 
efficiency  although  at  a  cost  greater  than  that  of  a  standard  flat  plate  system. 

Table  6  tests  sensitivity  to  the  effects  of  financing,  zero  interest  utility 
loans  and  fuel  price  escalation  rates.    The  test  case  system  is  not  economic  when 
financed  with  a  high  interest  loan.    This  benefit-cost  analysis  is  very  sensitive  to 
the  discount  rate  and  fuel  price  escalation  rates  chosen.    DNRC  explicitly  assumed 
that  fuel  price  escalation  rates  will  be  in  addition  to  inflation. 


Table  6.  -  Sensitivity  to  Financing  Method 
and  Fuel  Price  Escalation  Rate 
Natural  Gas  Electric 


Optimal 

ipi 

Net 

ft2 

•  pi 

Net 

ft2 

Savings 

Savings 

Reference  Case 

64 

.53 

$-242 

96 

.70 

$  835 

Financed 

53 

.46 

-477 

86 

.66 

349 

O-Interest  Loan 

86 

.64 

354 

107 

.74 

1,931 

1  percent  FP  escal.rate 

53 

.46 

-425 

86 

.66 

506 

5  percent  FP  escal.rate 

86 

.64 

229 

107 

.74 

1,618 

Reid  and  Hendricks  (1977)  noted  the  pronounced  effect  of  the  federal  tax  credit 
on  the  financial  attractiveness  of  solar.    This  provision  of  the  1978  Conservation 
Tax  Act  is  slated  to  expire  in  1984  and  its  renewal  is  uncertain.    Table  7  explores 
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effects  of  a  discontinued  federal  tax  subsidy  and  also  of  the  effect  of  increasing 
the  Montana  state  tax  credit  to  20  percent  of  initial  costs.    Results  of  solar 
competing  with  electricity  price  at  its  marginal  cost  $.07/kWh  (Duffield,  1979)  are 
also  shown.    This  measure  is  important  because  while  the  individual  choice  to  . 
install  a  solar  collector  is  based  on  the  average  costs  of  fuels,  the  actual  benefit 
to  society  is  that  the  SHW  system  displaces  the  need  of  new  electricity  generation 
capacity.    The  cost  of  this  new  capacity  -  marginal  cost  -  typically  is  much  higher 
than  the  average  cost. 


Table  7.  -  Sensitivity  to  Tax  Subsidy/Marginal  Cost  Pricing 
Natural  Gas  Electric 


Optimal 

ipi 

Net 

ft2 

Net 

ft2 

Savings 

Savings 

Reference  Case 

64 

.53 

$-241 

96 

.70 

$835 

No  Federal  Tax 

Credit 

21 

.21 

-1229 

53 

.47 

-566 

20  percent  State 

Tax  Credit 

86 

.64 

348 

107 

.74 

1546 

Marginal  Cost 

Electricity 

118 

.77 

4637 

Without  the  federal  tax  credit  solar  cannot  compete  economically  with 
conventional  energy  sources.    A  20  %  state  tax  credit  does  improve  the  financial 
attractiveness,  especially  for  conversion  from  natural  gas.    If  electricity  were 
priced  at  its  marginal  cost,  the  reference  case  system  would  provide  the  greatest 
present  value  savings  of  any  sensitivity  performed. 
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Uncertainty  Analysis 


In  the  economic  analysis  of  solar  options,  several  economic  parameters  are 
estimated  with  varying  degrees  of  uncertainty.    These  include  the  discount  rate, 
inflation  rate  and  fuel  price  escalation  rates.    Future  energy  price  projections  are 
especially  uncertain,  yet  they  have  profound  effect  on  the  results  of  the  cost 
benefit  analysis.    The  problems  with  forecasting  future  energy  prices  are 
demonstrated  in  the  following  examples. 

In  a  1976  report  from  the  Montana  Energy  Advisory  Council,  Stroup  (1976) 
projected  Western  Montana  natural  gas  prices  through  the  year  1990.    At  that  time 
(end  of  1976)  the  residential  price  of  natural  gas  was  $1.66  per  MCF;  the  report 
projected  the  1980  gas  price  at  $2.34  in  constant  1976  dollars.    At  the  close  of 
1980  the  price  of  natural  gas  was  $4.07/MCF  -  $3.03  when  deflated  to  1976  constant 
dollars.    Projections  were  made  that  the  price  of  gas  would  increase  at  an  annual 
rate  of  8.6  percent  per  year  -  certainly  not  a  conservative  estimate.    The  actual 
annual  percent  growth  rate  was  15  percent!    Furthermore,  Stroup  assumed  that  growth 
in  natural  gas  prices  would  be  essentially  flat  in  the  1980s,  increasing  to 
$2.40/MCF  (1976  dollars)  by  1990.    This  point  actually  was  passed  in  1979. 

In  a  companion  report  Duffield  (1977)  predicted  that  average  electricity  prices 
would  increase  at  about  2  percent  in  real  terms  over  the  period  1973-1990. 
Specifically,  in  the  base  case,  a  15.2  percent  cumulative  increase  in  constant  1976 
dollars  was  forecast  for  1973-1980  period  corresponding  to  a  2.36  percent  annual 
percent  growth  rate.    Converting  residential  electric  rates  into  1976  dollars,  the 
actual  commulative  rate  of  increase  was  11.8  percent,  a  1.80  percent  annual  growth 
rate.    This  is  good  agreement  with  actual  events,  but  does  overestimate  the  real 
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^  .   In  either  case  the  predictions  were  in  error.    Stroup  correctly  anticipated 
large  increases  in  natural  gas  prices,  but  not  nearly  large  enough.    On  the  other 
hand,  Duffield  made  a  conservative  guess  for  electricity  prices  but  not  conservative 
enough.    The  probability  of  forecasting  future  energy  prices  with  accuracy  seems 
low. 

Analysis  for  a  period  of  twenty  years  is  even  more  uncertain  than  the  two 
short-term  cases  just  examined.    The  DNRC  model  was  a  long-term  real  fuel  escalation 
rate  of  2.5  percent.    A  short-term  (ten  year)  escalation  rate  of  4.0%  for  natural 
gas  and  6.0%  for  electricity  is  assumed.    For  purposes  of  a  sensitivity  analysis,  we 
could  assume  that  the  fuel  price  escalation  rate  is  uncertain  such  that  the  actual 
value  for  »e'  over  the  period  may  vary  ±  2  percent  from  the  value  DNRC  used. 

* 

We  will  also  assume  that  the  individual  discount  rate  is  uncertain  such  that  the 

actual  value  of  "d"  over  the  period  may  vary  ±  1%  from  the  value  used  (2.78%) .  The 

uncertainty  in  either  variable  is  approximated  by: 

A  NS  =  _3_NS_  AXi 
9  Xi 

Where:      A  NS  =  change  in  the  present  value  of  net  savings 

9  NS  =  differential  of  net  savings  over  the  period 

9  Xi  =  differential  value  of  the  parameter  (discount  rate  or  fuel  price 

escalation  rate) 

A  x^=  the  change  in  "x^"  examined 

For  the  analysis,  the  same  reference  case  system  (four  person  household,  in 
Helena  with  an  80  square  foot  collector  array)  is  used.    Electricity  savings  in  the 
first  year  amounts  to  $157.    Use  of  the  reference  case  economic  values  yields  a  net 
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savings  of  $804  over  the  period.    Uncertainty  in  the  discount  and  escalation  rates 
affect  the  value  of  "UPW*"  as  shown  in  the  equation  on  page  11.    As  a  result  the 
uncertainty  in  net  savings  generated  by  the  discount  rate  is  such  that  the  net 
savings  may  vary  ±  $333.    The  same  calculation  for  the  fuel  price  escalation  rate 
yields  a  ±  $665  uncertainty.    The  maximum  uncertainty  generated  by  variations  in  the 
two  parameters  is  simply  the  sum  of  the  two:  +$998. 

However,  based  on  Brandemuehl  and  Beckman  (1979),  the  "most  probable" 
uncertainty  in  net  savings  is  given  by: 


so  that  the  most  probable  uncertainty  generated  by  the  two  parameters  =  ±  $744. 
Since  the  base  case  net  savings  is  only  $804,  the  analysis  indicates  that  with  only 
a  two  percent  uncertainty  in  future  fuel  prices  and  one  percent  uncertainty  in  the 
discount  rate,  the  actual  life  cycle  savings  may  range  from  $1,548  to  only  $61. 
This  ranges  from  excellent  savings  to  almost  nonel    This  underscores  the  sensitivity 
of  the  long-term  savings  predicted  by  any  SHW  system  simulation  to  even  small 
variations  in  exceedingly  uncertain  future  fuel  prices  and  discount  rates.  This 
uncertainty  is  true  not  only  for  solar  investments  but  also  for  those  involving 
conventional  fuels. 

Ranking  Solar  Thermal  of  Thirty  Montana  Sites 

Fcwlkes  (1981)  has  measured  insolation  and  climatic  data  for  30  Montana 
locations.  SHW  analysis  was  performed  for  each  of  these  sites  to  establish 
geographic  variation  in  solar  potential.    The  base  case  assumptions  (four  person 
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Ranking  Solar  Thermal  of  Thirty  Montana  Sites 

Fowlkes  (1981)  has  measured  insolation  and  climatic  data  for  30  Montana 
locations.    SHW  analysis  was  performed  for  each  of  these  sites  to  establish 
geographic  variation  in  solar  potential.    The  base  case  assumptions  (four  person 
household,  average  hot  water  use)  were  used  for  the  ranking  procedures.    The  ranking 
results,  based  on  delivered  solar  energy,  are  presented  in  Table  8  and  Figure  2.  An 
80-square-foot  collector  array  is  assumed.    Although  the  model's  absolute  accuracy 
is  uncertain,  the  ordinal  ranking  is  at  least  approximately  accurate.    The  major 
caveat  in  this  assessment  is  the  short-term  nature  of  the  Fowlkes  data.  Regardless, 
the  measured  Fowlkes  data  is  more  believable  than  SQLMET  estimates  previously 
available  from  the  U.S.  Weather  Service. 
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Table  8. 

Solar  Thermal  Ranking  of  30  Montana  Sites 


Solar  Heat 


"Fractionn " 

Collected^ 

Rank 

.71 

16.1 

1 

pari  nwi~nn 

.68 

15.5 

2 

ij  t-W  J.  O  L.       1  i 

68 

15.5 

3 

JUIQaJl 

•  DO 

15  5 

A 

Oil  1  i  rinc 

.  D 

15  4 

5 

Mil  c»g  Pn  +-\7 

68 

15.4 

6 

•  DO 

K  A 
1J  •  ** 

7 

fi7 

•  D  / 

n 

fi7 

•  D  / 

IS  4 

J-  -J  •  •* 

g 

biabyuW 

K7 
•  0  / 

J. -J  •  z 

l  n 

unoueau 

•  DO 

X3«U 

1 1 

vjreat  rails 

•  DD 

10  .u 

X4 

fx  en tywooa 

.03 

10. U 

13 

west  leiiowsuone 

.bo 

1A  R 
14  .  o 

~\  A 
14 

UiIUlXO 

•  DO 

1  K 

U  iv  xny  bt.  un 

fit; 

.00 

1  A  7 
14.  / 

XO 

.DD 

14.0 

X  / 

KA 
•  04 

J.4  •  D 

1  R 

XO 

11  CiV  J.  c 

14  4 

19 

DUUwC 

14  4 

20 

Del  CI  la 

•  UJ 

111  J 

21 

•  UJ 

14  1 

22 

Fort  Benton 

.63 

14.2 

23 

Anaconda 

.61 

13.9 

24 

Hamilton 

.60 

13.6 

25 

Kalispell 

.56 

12.9 

26 

Poison 

.56 

12.7 

27 

Missoula 

.56 

12.7 

28 

Thompson  Falls 

.54 

12.1 

29 

Libby 

.53 

12.1 

30 

Notes:       80  ft2  solar  collector;  FRUL  =  .68,  FR  xa=  .72 

4  person  occupancy,  average  hot  water  use,  tank  set  point  of  140°F 

1  percentage  of  thermal  load  provided  by  solar  collectors 

2  106  Btus/year 
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Figure  2 

Montana  Solar  Hot  Water  Thermal  Performance  Ranking 
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The  Economic  Ranking 


The  economic  ranking  of  the  thirty  Montana  solar  sites  varies  appreciably  from 
the  thermal  potential  ranking.    This  is  because  current  prices  of  energy  differ 
across  the  state  depending  on  the  utility  or  local  cooperative  supplier. 
Electricity  prices  have  the  greatest  range  of  variance.    Hydroelectric  power  in  the 
Flathead  Irrigation  Project  service  area  is  only  $.019/kWh  while  thermally  produced 
electricity  in  the  MDU  service  area  is  $.069/kWh.    Natural  gas  is  not  available  in 
some  sites. 

The  economic  ranking  is  presented  in  Table  9  and  Figure  3.    There  is  greater 
uncertainty  in  this  analysis  than  the  physical  one.    Only  low  confidence  is  possible 
in  estimates  of  the  real  energy  price  escalation  rates  and  how  they  might  vary  from 
one  utility  to  another. 
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Table  9 

Solar  Economic  Ranking  of  30  Montana  Sites 
(Arranged  in  order  or  rank) 


ELECTRICITY  NATURAL  GAS 


LilSt.  or 

Net 

List  of 

Net 

Towns 

Savings 

Rank 

Towns 

Savings 

Ban] 

niies  city 

$4,396 

1 

Great  Falls2 

$  104 

1 

Glasgow 

4,344 

2 

Dillon 

24 

2 

rientywooa 

4,233 

3 

Jordan 

-  60 

3 

Gienaive 

4,044 

4 

Lewi st own 

-  63 

4 

Sidney 

3,931 

5 

Harlowtown 

-  72 

5 

Kalispelli 

1,412 

6 

Red  Lodge 

-  79 

6 

Dillon 

1,226 

7 

Browning 

-  80 

7 

LlDDy^ 

1,156 

8 

Billings 

-  83 

8 

Jordan 

1,103 

9 

Colstrip 

-  86 

9 

bewistown 

T  AAA 

1,099 

10 

Choteau 

-149 

10 

naricwtown 

1,086 

11 

Livingston 

-192 

11 

Rea  Lodge 

1,079 

12 

Ennis 

-196 

12 

Billings 

1,071 

13 

Bozeman 

-223 

13 

loi strip 

1,065 

14 

Butte 

-232 

14 

W.  Yellowstone 

1,040 

15 

Havre 

-242 

15 

iwjioteau 

975 

16 

Helena 

-267 

16 

crear  tans 

964 

17 

Fort  Benton 

-282 

17 

.Livingston 

ft!  1 

911 

18 

Anaconda 

-330 

18 

rjinlS 

908 

19 

Hamilton 

-371 

19 

Bozeman 

867 

20 

Miles  Citv 

-47  ^ 

HID 

on 

Butte 

855 

21 

Kalispell 

-488 

21 

Havre 

841 

22 

Glasgow 

-490 

22 

Helena 

804 

23 

Missoula 

-517 

23 

Fort  Benton 

780 

24 

Plentywood 

-524 

24 

Anaconda 

714 

25 

Glendive 

-581 

25 

Hamilton 

652 

26 

Thompson  Falls 

-615 

26 

Missoula 

439 

27 

Sidney 

-617 

27 

Thompson  Falls 

296 

28 

Libby 

N/A 

Browning 

264 

29 

W.  Yellowstone 

N/A 

Poison 

-  983 

30 

Poison 

N/A 

Notes:    80  ft2  solar  collector;  FrO"l=.68,  FRtq  =.72 

4  person  occupancy,  average  hot  water  use,  energy  cost  summary  from  Table  3. 

1-  mean  of  rate  tiers.    See  Table  2. 

2-  Great  Falls  gas  price,  all  others  are  MPC,  MDU. 
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Figure  3 

Montana  Solar  Hot  Water  Economic  Performance  Ranking 
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Summary 

Based  on  the  simulation  model  assumptions,  solar  hot  water  heating  is 
competitive  with  electricity  in  Montana  but  is  not  competitive  with  natural  gas  hot 
water  heating.    Optimal  collector  size  is  about  64-96  square  feet  for  a  family  of 
four  -  the  larger  figure  for  solar  against  electric  heating.    The  greatest  problem 
is  the  long  payback  time  involved  because  of  the  large  initial  costs. 

Factors  such  as  larger  households  or  other  influences  that  increase  hot  water 
use  tend  to  favor  solar  investments.    Self-built  systems  are  considerably  cheaper 
than  the  commercial  variety  although  resulting  system  efficiency  may  be  lower. 
Owner  installed  systems  seem  to  offer  no  significant  advantage  given  the  risk  of  a 
poor  installation.    (See  Appendix  A) 

Solar  hot  water  heating  is  not  economically  viable  in  Montana  without  the 
federal  tax  credit.    An  increased  state  tax  credit  would  increase  solar's  economic 
competitiveness  -  especially  for  natural  gas  users. 

Any  economic  evaluation  model  is  particularly  susceptible  to  distortions 
stemming  from  selection  of  the  discount  rate,  fuel  price  escalation  rates  and 
assumptions  regarding  future  inflation.    Special  care  should  be  taken  in 
specification  of  these  values.    The  effect  of  the  uncertainty  of  these  values  has 
been  illustrated.    A  2  %  variation  in  the  fuel  price  escalation  rate  and  a  1  % 
uncertainty  in  the  discount  rate  could  have  the  effect  of  negating  any  solar 
economic  benefit  or  conversely,  would  triple  life  cycle  solar  savings. 
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Lower  borrowing  rates  and  higher  inflation  favors  solar  investments  in  that 
effective  discount  rates  are  reduced.    Study  of  the  impact  of  various  policy 
oriented  incentive  programs  for  solar  investment  is  easily  accoumpl ished  with  the 
model.    If  electricity  were  priced  at  its  marginal  cost  ($.07/kWh) ,  solar  would 
currently  be  very  competitive  in  all  locations. 

Fowlkes  (1982) ,  having  measured  the  variation  in  solar  flux  at  30  sites  in 
Montana,  found  that  eastern  sites  receive  greater  amounts.    The  economic 
attractiveness  of  solar  hot  water  heating  in  Montana  tended  to  follow  these 
findings.    The  best  sites  are  east  of  the  Continental  Divide.    The  worst  are  located 
in  western  Montana.    A  typical  80  square  foot  collector  array  would  provide  53  to  71 
percent  of  the  annual  hot  water  needs  in  the  various  locations.    First  year  savings 
ranged  widely  from  about  $70  in  Poison  to  over  $300  in  Miles  City. 

The  economic  ranking  of  the  sites  differs  considerably  from  the  thermal  ranking 
because  of  the  large  effect  of  local  conventional  energy  prices.    The  higher  natural 
gas  rates  in  the  Montana  Power  Company  system  tended  to  favor  those  sites  for  sclar 
over  sites  served  by  Montana-Dakota  Utilities.    However,  Montana-Dakota  Utility 
electric  rates  are  currently  the  highest  in  Montana  at  6.9  cents  per  kwh.    The  solar 
sites  in  this  area  (Glasgow,  Plentywood,  Miles  City)  had  the  shortest  payback 
periods  examined  (7  years) .    The  inverted  rate  structure  of  the  PP&L  area 
(Kalispell,  Libby)  has  the  effect  of  greatly  favoring  solar  as  monthly  electrical 
consumption  increases  above  1,000  kwh,  as  is  the  case  with  homes  using  electric 
resistive  heating.    The  attractiveness  of  SHW  in  the  co-ops  examined  varied 
depending  on  their  current  rate  structure. 
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The  greatest  economic  barrier  to  utilization  of  the  sun's  energy  is  the  high 
initial  cost  of  the  systems  and  the  comparatively  low  avoided  cost  of  conventional 
energy  consumption.    One  may  assume  that  the  full  potential  for  reducing  the  initial 
costs  of  solar  systems  has  not  yet  been  realized.    For  instance,  development  of  new 
polymer  film  solar  collectors  promises  to  reduce  materials  and  manufacturing  costs 
of  collectors  by  an  order  of  magnitude  (Wilhelm  and  Andrews,  1982) .    Advertised  kit 
systems  seem  to  offer  substantially  lower  prices.    Advanced  evacuated  tube  systems 
promise  high  performance,  although  costs  must  be  reduced  somewhat  to  be  commercially 
attractive.    Given  this,  we  may  expect  the  initial  cost  of  the  systems  to  decrease 
in  real  terms. 

There  is  evidence  that  consumer  time  horizons  for  energy  investments  are  not 
appreciably  greater  than  five  years  (Cunningham  and  Lopreato,  1977) .    As  we  have 
seen,  even  the  most  favorable  analysis  has  larger  paybacks  than  this.    However,  most 
large  scale  hydro-power  sites,  the  lowest  cost  energy  sources  in  Montana,  have  been 
exploited.    New  additions  to  generation  capacity  will  usually  consist  of  thermal 
conversion  facilities  with  higher  priced  electricity.    Further,  the  real  price  of 
conventional  fuels  is  likely  to  increase  over  the  long  term,  and  coupled  with  any 
reduction  in  manufacturing  costs,  this  would  greatly  increase  the  life-cycle  savings 
of  solar. 

Solar  hot  water  heating  is  now  only  provisionally  economic  in  Montana.    It  is 
probable  that  long-term  real  fuel  price  increases  will  bring  the  technology  into 
common  use  in  the  state  by  the  twenty-first  century. 
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APPENDIX  A 
TOE  FCHART  ALGORITHM 


Two  dimensionless  parameters  are  calculated  -  "X"  and  "Y"  that,  when  correlated, 
yield  "F"  -  the  monthly  solar  fraction: 

X  =  h°E^JT^f-  fk)  At 
L 

Y  -  AoEr(  fa  )HtS 


L 

where.       Ac  _  co^e^or  area 

FR  =  collector  -  heat  exchanger  efficiency  factor 

UL  =  collector  heat  loss  coefficient  (W/M2C°) 

At  =  number  of  seconds  in  the  month 
T 

a  =  monthly  average  ambient  temperature  (C°) 

Tref  =  reference  temperature  (100°C) 

L  =  monthly  heating  load  (Joules) 

Ht  =  monthly  average  radiation  incident  on  the  collector  surface 
(J/m2) 

N  =  days  per  month 

•f^x  -  monthly  average  transmittance-absorbance  product 


The  products  FrUL  and  FrTa  are  determined  from  controlled  collector  tests 
and  then  modified  by  a  heat  exchanger  penalty  factor:  f 

I  F' 

which  is  assumed  to  be  .96  for  typical  units.    Tau-alpha  is  modified  by  the  collector 
incident  angle  to  yield  a  corrected  absorbance  transmittance  product.    This  is 
assumed  to  be  .94  for  standard  double  glazed  designs. 
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To  use  FCHART  for  hot  water  heating  systems,  an  additional  correction  factor  for 
"X"  must  be  calculated.  The  main  water  supply  temperature  and  minimum  acceptable  hot 
water  temperature  both  effect  collector  thermal  losses.  is  calculated: 

Xc  =  (11.6  ±  1.18TW  ±  3.86  Tm-^LJ^) 
X  (100  -  Ta) 

Where:  =  multiplier  for  Xc 

X 

Tw  =  minimum  acceptable  hot  water  temperature  (C°) 

Tm  =  mains  water  supply  temperature 

Tc  =  monthly  average  ambient  temperature  (C°) 

A  correction  factor  must  also  be  calculated  for  the  monthly  load.    This  factor 

must  take  into  account  thermal  losses  from  the  auxiliary  tank,  which  are  expressed 

■ 

as: 

Laux  -  UA  (Tw  -  T'a) 

Where:       U  =    thermal  loss  coefficient  (W/m2) 
A  =   area  of  tank 
Tw  =  temperature  of  hot  water 

T*a  =  ambient  air  temperature  around  storage  tank 

In  the  model  the  auxiliary  tank  is  assumed  to  be  of  standard  design  with  R-10 
external  insulation:  "u"  which  accounts  for  both  tank  and  plumbing  losses  is  given  to 
be  .62  W/m2  (electric)  and  .82  w/m2  (gas),  A  is  3.21  m2,  Tw  =  60°C  and  T'a  is  16°. 
The  corrected  load  factor  is: 

Lc,i  =  Li  +  Laux 
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The  results  are  two  dimensionless  variables  Xc  and  Y  that  show  the  reference 
collector  heat  loss  and  solar  gain  over  the  load.    Based  on  these  values,  the  monthly 
solar  fraction  is  derived  mathematically  using  a  power  series  correlation  for  liquid 
systems: 

Fmonth  =  L029Y  -  .065X  -  .245Y2  +  .0018X2  +  .0251Y3 

■ 

The  fraction  of  the  annual  load  supplied  by  solar  is  the  sum  of  monthly  solar 
contributions  ( i  FjL-^) ,  divided  by  the  annual  load  (  £  L^) . 

12 

Fa  =      I  _£iLi 

■     Li  I 
i=l 

The  fraction  may  be  examined  for  different  collector  areas,  locations,  and 
collector  types.    These  configurations  may  then  be  compared  in  the  life-cycle  costing 
methodology. 

The  FCRART  algorithm  has  demonstrated  less  predictive  capability  when  applied  to 
non-conventional  collectors,  improperly  installed  units  and  homebuilt  systems  of 
unknown  specification  (Winn,  1981) .    Winn  found  that  FCRART  typically  overestimates 
the  fraction  of  these  systems  by  22  percent.  Consequently,  the  monthly  fraction 
predicted  for  homebuilt  systems  has  been  uniformly  reduced  by  20  percent  for  the 
purposes  of  the  analysis.    Owner  installed  systems  have  been  decreased  by  10  percent 
because  of  the  risk  of  poor  installation. 
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Appendix  B 

1)  The  size  of  the  demand  for  hot  water  in  a  household  of  four  is  calculated  as 
follows: 

■ 

Based  on  data  collected  by  Bavir  et  al.  (1981)  per  person  water  consumption 
is  best  described  by  an  exponential  decay  function  according  to  the 
equation: 

Q  =  P*  (A  *  P"x) 
Q  =  P*  (23.4  *  (P)".18 
(R2  =  .725) 

Where: 

Q  =  gallons  of  hot  water  used  per  day 
P  =  number  of  people  in  the  household 
x  =  fitted  exponential  coefficient 
A  =  the  regression  intercept 

2)  Cold  water  line  temperature  varies  from  40°  F.  in  late  winter  to  60°  in  early 
fall.    This  is  based  on  a  year-long  measurement  of  mains  water  temperature  by 
month  in  two  western  Montana  sites.    A  Weksler  Fahrenheit  thermometer  was  used  to 
take  the  measurements  in  Missoula  and  Helena  on  the  15th  of  each  month.  The 
thermometers  were  calibrated  monthly.    The  relationship  is  a  lagged-linear 
expression  assuming  the  residual  values  to  be  serially  independent. 

The  function  is  of  the  form: 
%  -  .38  1^!  +33.4 
(R2  =  .991) 

Where:       TA_^  =  ambient  temperature  in  previous  month 
Tw  =  water  temperature  in  current  month 
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3)  The  hot  water  tank  set  temperature  is  140°F  (Palmiter  and  Baylin,  1982) . 

4)  The  average  heat  necessary  to  heat  this  quantity  of  water  is  calculated: 

Li  =  Q  *  VT  (1h  -  Tc) 
Where:       L  =  daily  thermal  load  (Btus/day)  in  month  "i" 
Q  =  gallons  consumed  per  day 
VT  =  volumetric  weight  of  water  (8.33  lbs/gal.) 
Th  =  hot  water  set  temperature 
Tc  =  mains  water  temperature 

Monthly  fuel  demand  is  calculated: 
G  =  jgj  *  b) 
(HC    *  EC) 

G  =  monthly  fuel  consumption  (mcf ,  kwh) 
Where:      D  =  number  of  days  in  month  "i" 

HC  =  heat  content  of  natural  gas  at  an  altitude  of  4,000  ft.  (924,000 
Btu/ft3)  or  1  kWh  of  electricity  (3413  Btus) . 

EG  =  fractional  efficiency  in  delivering  heat  content  of  natural  gas  or 
electricity  to  water  including  standby  losses,    (natural  gas:  .47 
no  external  insulation,  .55  insulated  to  R-ll;  electricity:  .76  no 
external  insulation,  .85  insulated  to  R-ll  Farahan,  1977,  Ek,  1981). 
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380 

390 

400 

410 

420 

430 

450 

470 

480 

500 

501 

502 

503 

504 

506 

510 

51 1 

512 


REM  * 
REM  * 
REM  * 
REM  * 


SOLAR  HOT  WATER 
ANALYSIS  MODEL 


REM  ************************ 


REM 
REM 
TEXT 
JA$  = 
MR*  = 
MY*  m 
JL*  = 
SP*  = 
NV*  ■ 
D*  = 
DEF 


"JANUARY" : FB*  =  "FEBRUARY" 
" MARCH ":AP*  =  "APRIL" 
"MAY":JE*  ■  "JUNE" 
"JULY" : AG*  =  "AUGUST" 
"SEPTEMBER" : OT*   =  "OCTOBER" 
"NOVEMBER": DC*  =  "DECEMBER" 
CHR*    (4)  :    REM     CTRL  D 

FN  C(X)    ■    (X   -  32)    /    1.8:    DEF     FN  F(X)    =   X   *    1    8  +  32 
E>  I M  F I  $  ( 30 )  , A ( 40 )  ,  N  ( 1 00 )  ,  V  ( 1 00 )  ,  P  ( 1 00 ) 

HOME   :    INVERSE   :    GOSUB  21380:    PRINT   "  *   D0  YOU  HAVE  A  PR TNT 

ER?";:    GOSUB  21380:    GOSUB  21380:   NORMAL   :  PRINT 
PRINT   -        1)    YES":    PRINT   "       2)   NO":  PRINT 
CV  =  2s PF  -  OsLO  -  Is  HI  =  2:   GOSUB  39300 
IF  CV  =   1   THEN  PF  * 1 

IF  PF  THEN     HOME   :    INVERSE   :   GOSUB  21380:   PRINT   "   *IN  WHICH  qi  nr 
IS  THE   INTERFACE  CARD-;:    GOSUB  21380:    GOSUB  2 1 380 :    NORMAL  f  PRINT 

IF  PF  THEN  CV  =   1:L0  =   Is  HI   -  8:    GOSUB  39300: PS  =  CV 

GOSUB   10400:    REM   I NT I T I AL I ZE  CITY  FILE  NAMES 

GOSUB   10050:    REM     READ  CITY  FILE 

GOSUB   10530:    REM     READ  DEFAULTS 

REM 

REM  *  REASSIGN  VARIABLES  * 
REM 

Al  =  A(l) : A2  =  A (2) 
A3  =  A  (3)  :  A4  m  A  (4) 
A5  A  (5)  :  A6  =  A  (6) 
A7  =  A(7) : A8  =  A<8) 
A9  ■  A (9) : AO  =  A (10) 


AY  =  A(ll) 
A (10)  =  AO 
II  =  A(13) 
13  m  A ( 15) 
15  ■  A(17) 
17  =  A(19) 
19  =  A(21) 
IA  «  A (23) 
IF   IA  +  IB 


II  -  II  * 
15  =  15  * 
19  =  19  * 
GOTO  520 
II  =  II  * 
15  =  15  * 
19  =    19  * 


AZ  =  A ( 12) 
A (11)  =  AY 
12  =  A(14) 
14  =  A(16) 
16  =  AU8) 
18  =  A (20) 
10  =  A (22) 
IB  =  A(24) 

1550  THEN 


.94:12  =   12   *  , 
1.13:16  =    16  * 
1 .  04 :  1 0  =    1 0  * 

. 98: 12  ■    12   *  , 
1  .  12: 16  =   16  * 
1 . 04 : 10  m    10  * 


510 

.96:13  =   13   *    1.02:14  =  14   *  1.08 

1.15:17  =   17   *    1.14:18  =   IB   *  1.1 

.98:IA  =   IA   *    .95:IB  =  IB   *  .93 

.99:13  =   13   *    1.04:14  =  14   *  1.08 

1.15:17  -   17   *    1.14:18  =   18   *  1.1 

.99: I A  =    I A   *    .96:  IB  =  IB   *  .96 


520 

11   =   11  * 

. 01 13538: 12  = 

12 

*  . 

<_>  1  1  •_*'  wJ  •  J"  O 

530 

13  =   13  * 

.011 3538 : I 4  = 

14 

(.)  1  1  OJOO 

540 

15  =   15  * 

. 01 13538: 16  = 

16 

(.)  1  1  OJOD 

550 

17  =  17  * 

. 0 1 1 3538 : I 8  = 

18 

V  1  lJJ O O 

560 

19  =   19  * 

. 01 13538: 10  = 

10 

(.>  1  1  OJOO 

570 

IA  =   IA  * 

. 01 13538: IB  = 

I  b 

* 

t~\  1  1  IR^R 
Ul 1 OJOO 

620 

Al   -  <A1 

-  32)    /  1.8:A2 

= 

<  A2 

—     32)  1 

=    <A4  - 

32)    /  1-8 

630 

A5  =  <A5 

-  32)    /  1.8:A6 

( A6 

—      32)  / 

=    <A8  - 

32)    /  1.8 

640 

A9  =  (A9 

-  32)    /    1.8: AO 

—     O  i.  1  / 

=    (AZ  - 

32)    /  1.8 

650 

HOME 

660 

REM 

670 

REM 

680 

TEXT   :    HOME   :  INVERSE 

:    GOSUB  21380: 

BER  OF  PEOPLE";:  GOSUB 

21 

380 

:  PRINT 

1.8: A3  =  (A3  -  32)  /  1.8:A4 
1.8:A7  =  (A7  -  32)  /  1.8:A8 
1.8:AY  -   < AY  -  32)    /  1.8:AZ 


PRINT 


GOSUB  21380:    GOSUB  21380:    NORMAL   :  PRINT 


*  ENTER  THE  NUM 
IN  YOUR  HOUSEHOLD"; 
PRINT- :  PRINT 


690 
700 
710 
720 


1:HI   -    100:    GOSUB  39300: ZP  =  CV:P  =  ZP 


730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 

860 
870 
880 
890 


CV  =   ZP:L0  = 
REM 
REM 

HOME   :    INVERSE   :   GOSUB  21380: 
ERIZE";:   GOSUB  21380:   PRINT  ' 
21380:   GOSUB  21380:    NORMAL  : 
NORMAL 

MISERLY" 
AVERAGE" 
LAVISH":  PRINT 
=  liHX   =  3:  GOSUB 
THEN  Q  =    16.4   *  P 
THEN  Q  =  23.4   *  P 
THEN  D  -  30.4   *  P 


PRINT   "  *   HOW  WOULD  YOU  CHARACT 

YOUR  HOT  WATER  USAGE?";:  GOSUB 
PRINT   :    PRINT   :  PRINT 


PRINT  "1) 
PRINT  "2) 
PRINT  "3) 
CV  =  PY.-LO 
IF  PY  =  1 
IF  PY  =  2 
IF  PY  =  3 
EO  =  8.33 
REM 
REM 

INVERSE 

HOME   :  INVERSE 
MP" ; :    GOSUB  21 
CV  » 
REM 
REM 
HOME  : 
LING"  ; 
21380: 
900  PRINT 
910  PRINT 
920  PRINT 
930  CV   =  FA 
940  REM 
950  REM 
960     HOME  : 
ANCE" ; 
21380: 
PRINT  " 
PRINT  " 
PRINT  " 
CV  =  AA 


39300: PY  =  CV 
-  .18 


18 
18 


:    GOSUB  21380:    PRINT  " 
580:    GOSUB  21380:  NORMAL 


*  ENTER  HOT  WATER  TE 
PRINT   :    PRINT   :  PRINT 


FN  F(TH):LO  =  40: HI   =  212::    GOSUB  39300: TH  =     FN  C(CV) 


INVERSE   :    GOSUB  21380:    PRINT    "  *  HOW  ARE  YOU 

GOSUB  21380:    PRINT    "  THE  SOLAR  SYSTEM?";: 

GOSUB  21380:    NORMAL   :    PRINT   :    PRINT   :   PRINT  : 
1)    BUILDING/INSTALLING  YOURSELF" 
PRINT    "2)    BUYING  IT/COMMERCIAL 
PRINT   "3)    BUYING   I T/ I NSTALL I NG 
LO  =    1:HI    =  3:    GOSUB  39300: FA  = 


INSTAL 
GOSUB 


INSTALLATION" 
YOURSELF":  PRINT 
CV 


970 
980 
990 
1 000 
1 0 1 0 
1 020 
1030 


INVERSE   :    GOSUB  21380:    PRINT  " 

GOSUB  21380:  PRINT  "  THE 
GOSUB  21380:    NORMAL    :    FRINT    :  PRINT 

1)  A  BANK  OR  OTHER  LOAN" 

2)  USE    IT   AS   AN   ALTERNATIVE  INVESTMENT' 

3)  LOW   INTEREST  UTILITY  LOAN":  PRINT 
LO  =    1:HI    =  3:    GOSUB  39300: AA  =  CV 


*   HOW   WILL   YOU  FIN 
INSTALLATION?";:  GOSUB 
:  PRINT 


IF   AA  =    1  THEN 
IF  AA  =  2  THEN 
ZD  =  AA 


D  =  .046 
D  =  .0278 


1040 

IF  AA  = 

3 

THEN 

1 050 

Tl 

= 

.  38 

* 

AZ 

+ 

1060 

T2 

= 

.  38 

* 

Al 

+ 

1070 

T3 

= 

.  38 

* 

A2 

+ 

1 080 

74 

= 

.  38 

* 

A3 

* 

1 090 

T5 

= 

.  38 

* 

A4 

+ 

1 1 00 

T6 

.  38 

* 

A5 

+ 

1110 

T7 

= 

.  38 

* 

A6 

+ 

1120 

T8 

= 

.  38 

* 

A7 

♦ 

1130 

T9 

= 

.  38 

* 

A8 

+ 

1 1 40 

TO 

.38 

* 

A9 

+ 

1 1 50 

TA 

.  38 

* 

AO 

1 1 60 

TB 

.  38 

* 

AY 

+ 

1 165 

HOME 

:  INVERSE 

EC. 

PER  KWH " ; : 

c 

D  =  . 000000 1 
7.46 
7.46 
7.46 
7.46 
7.46 
7.  46 
7.46 
7.  46 
7.46 
7.46 
7.46 
7.46 

:    GOSUB  21380:    PRINT    " * ENTER  CURRENT  PRttf  nr  eh 
:.rPER  KWH";:    G0SUB  2.380=    GOSUB  2.380.    NORMAL  ^  PRINT  ^PR^NT 

CG!L°  "   -0l!H1SH  -   -12=    GOSUB  3<?300:CG  «  CV 

11/0  REM 
1180  REM 

1190     HOME   :    INVERSE   :    GOSUB  21380:    PRINT   »  *   HOW  DO  YOU  THINK  TM 

PRICE  OF";.-    GOSUB  21380:    PRINT   "  ELECTRICITY  ML^i™ » 

:    GOSUB  21380:    PRINT    -  OVER  THE  NEXT  20  YEAR'S^" ^-    GOBUB  51 

380:    GOSUB  21380:    NORMAL   :    PRINT   :    PRINT  PRINT 
1200     PRINT   "1)    SLOWLY  (17./YR)" 
1210     PRINT   "2)    MODERATELY  (37./YR)" 
1220     PRINT   "3)    RAPIDLY    <5'/./YR)  " 
1230     PRINT   "4)    VERY  RAPIDLY    (10%/YR)":  PRINT 
\lt°  C^       QE:L0  =   1« HI   =  4:    GOSUB  39300:  QE  =  CV 
1250     IF  QE  =   1    THEN  ES  =  .01 
1260     IF  QE  =  2  THEN  ES  =  .025 
1270     IF  QE  =  3  THEN  ES  =  .05 
1280     IF  QE  =  4   THEN  ES  =  .10 

1282     HOME   :    INVERSE   :    GOSUB  21380:    PRINT   » *  ENTER  FED  TAX   CP  PPT  FrtR 

ANALYSIS" ; :    GOSUB  21380:    GOSUB  21380:    NORMAL   :    PRI^   t    PRINT   . PRINT 

1284  CV  =  TX;LO  =  0:HI    =  75:    GOSUB  39300: TX   =  CV 
1290  WT  =  Q   #   3.78544   *   P   *  4190 

1300  LI   =    (WT  *  (TH  -  Tl)  *   31)    /    in    •  9 

1310  L2  -    (WT   *  (TH  -  T2)  *   28)    /    10  9 

1320  L3  =    (WT   *  (TH  -   T3)  *   31)    /    10  •-  9 

1330  L4  =    (WT   *  (TH  -  T4)  *   30)    /    10  •  9 

1340  L5  ■    (WT   *  (TH  -  T5)  *   31)    /    10  9 

1350  L6  =    (WT   *  (TH   -   T6)  *   30)    /    l0     ■  9 

1360  L7  -    (WT   *  (TH  -  T7)  *  31)    /    10  9 

1370  L8  =    (WT  *  (TH  -  T8)  *  31)    /   10  9 

1380  L9  m   (WT  *  (TH  -  T9)  *  30)    /   10   '  9 

1390  LO  =    (WT   *  (TH  -  TO)  *   31)    /    ID   ••  9 

1400  LA  =    (WT  *  (TH  -   TA)  *   30)    /    10  '  9 

1410  LB  =    (WT   *  (TH  -  TB)  *   31)    /    10  '-  9 
1420  HC  ■  .0036 
1430  REM 
1440  REM 

1450     HOME    :    INVERSE    :    GOSUB  21380:    PRINT   "  *    IS  YOUR  HOT  WA 

rn^t  213801    PRINT    "  TANK   INSULATED-::  GOSUB 

21380:    GOSUB  21380:    NORMAL   :    PRINT    :    PRINT    :  PRINT 
1460     PRINT    "1)  YES" 
1470     PRINT   "2)    NO":  PRINT 

1480  CV  =   Z:LO  =   It HI    =  2:    GOSUB  39300: Z   =  CV 
1490      IF   Z   =   1   THEN  EG  =  .85 
1500      IF  Z   =  2  THEN  EG  =  .76 


1510 

EH  = 

HC 

EG 

1520 

G2  = 

L2 

/ 

EH:  Gl 

L 1 

/ 

fc.n 

1 530 

G3  = 

L3 

/ 

EH:  G4 

1  A 

L4 

/ 

tn 

1540 

G5  = 

L5 

/ 

EH:  G6 

L6 

/ 

1 550 

G7  = 

L7 

/ 

EH:  G8 

= 

L8 

/ 

cri_i 
tH 

1560 

G9  = 

L9 

/ 

EH:  GO 

/ 

trl 

1570 

GA  = 

LA 

/ 

tH:  bb 

ru 

1 5B0 

REM 

1 590 

REM 

1 600 

HOME  : 

INVERSE 

I 

GOBUB  2 

CTRICITY";:   GOBUB  213B0: 


1 605 

1610 

1 620 

1 630 

1640 

1 650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1745 

1750 

1 760 

1770 

1780 

1790 

1B00 

1810 

1 820 

1830 

1840 

1850 

I860 
1870 
1 880 
1890 


1 900 
1910 
1915 
1920 
1930 


1935 
1 940 
1 960 
1970 
1 980 


21380: 
LM  =  5: 
PRINT 
PR  I  NT 
PRINT 
PRINT 
PRINT 
PR  I  NT 
PRINT 
PRINT 
PRINT 
PRINT 
FRINT 
PRINT 
PRINT 
AG  =  Gl 
LM  -  0: 
PRINT 
PRINT 
UA 
ST 


GOSUB  21380:  NORMAL 
GOSUB  22120:  PRINT 
JA* 


1380:  PRINT 
PRINT  " 

PRINT 


*  MONTHLY  KWH  OF  ELE 
TO  HEAT  WATER";:  GOSUB 


TAB( 

20)  ; 

INT 

(Gl ) 

TAB  ( 

20 )  ; 

I  NT 

(G2) 

TAB  ( 

20)  ; 

INT 

(G3) 

TAB  ( 

20)  ; 

I  NT 

(G4) 

TAB  < 

20)  ; 

INT 

<G5) 

TAB  ( 

20)  ; 

INT 

(G6) 

TAB  ( 

20)  ; 

INT 

(G7) 

TAB  < 

20)  ; 

INT 

<G8) 

TAB  ( 

20)  ; 

INT 

(G9) 

TAB  ( 

20)  ; 

INT 

(GO) 

TAB  < 

20)  ; 

INT 

(GA) 

TAB  ( 

20)  ; 

INT 

(GB) 

FB*: 
MR*: 
AP*: 
MY*: 
JE*: 
JL*: 
AG*: 
SP*: 
OT*: 
NV*: 
DC* : 
:  PRINT 

+  G2  +  G3  +  G4  +  G5  +  G6  +  G7  +  G8  +  G9  +  GO  +  GA  +  GB 
GOSUB  22120 
PRINT   "ANNUAL  KWH  USE  =" ;    TAB (  19) 


LI 
L3 
L5 
L7 
L9 
LA 
ZL 
* 


PRINT 
.62  %  3  •  2 1 
(UA   *  30.4 
LI   +  ST:L2 


:  L3 

:  L5  ■ 

=  L7 

■■  L9  ■ 

'  LA 

=  (LI 
9479 


+  ST 


L4 
ST:  L6 
ST:  L8 
ST:LO 
ST:  LB 
L2  + 


HIT  ANY 

*  (TH  - 

*  24  * 

=  L2  + 
=  L4 
=  L6 
=  L8 
=  LO 
=  LB 
L3 


+ 
+ 
+ 


KEY  TO  CONTINUE 
16) 

3600)    /    10  ■'■  9 

ST 

ST 

ST 

ST 

ST 

ST 


I  NT 


(AG) 
GET  AK* 


+  L4  +  Li 


+  L6  +  L7  +  L8  +  L9  +  LO  +  LA  +  LB) 


GOSUB 

GOSUB 
P   -  ZP 

HOME  : 
HAVE"; : 
:  GOSUB 

PRINT 


5720 
5470 


INVERSE   :    GOSUB  21380:  PRINT 
GOSUB  21380:    PRINT  " 
21380:    GOSUB  21380: 
'1)    YES":    PRINT  "2) 


*  WOULD  YOU  LIKE  TO 
THE  SYSTEM  SIZED  FOR  YOU?": 
NORMAL   :    PRINT    :    PRINT    :  PRINT 
NO":  PRINT 


HI   =2:    GOSUB  39300: BQ  =  CV 


CV  =  BQ: LO  =  1 
XI   =   XI    *  5.67 
IF   BQ  =  2   THEN  2010 

HOME   :    INVERSE    :    GOSUB  21380:    PRINT  " 
HE  COLLECTOR" ; :    GOSUB  21380:    PRINT  " 
21380:    GOSUB  21380:  NORMAL 

GOSUB  22200: TM  =  Y5  +   is    PRINT   :  PRINT 
SA  -  1 

GOSUB  3870 
SA  =   SA  —  1 

TEXT   :    HOME    :    INVERSE   :    GOSUB  21380:  PRINT 
OR  AREA,    SQ  FT=";    INT    ( SA   *    10.7639);:  GOSUB 

NORMAL    :    PRINT    :    PRINT   :  PRINT 


PLEASE 
SIZE 

PR  I  NT 


STAND  BY  WHILE  T 

IS  OPTIMIZED"; :  GOSUB 


GOSUB  22100 


1  *  OPTIMAL  COLLECT 
21380:    GOSUB  21380: 


1999  SB  =  SA 

2000  PRINT  :  PRINT  :  PRINT  "HIT  ANY  KEY  TO  CONTINUE  .  .  .»,,  GET  AK< 
2010     HOUE   :    INVERSE   :    GOSUB  21380:    PRINT   »  *  WHAT  COL LFrrnp  fiDr. 

:U«60S^c21380!    PRINT   "  WOULD  YOU  uJeTSLSSSJ^bSS 

21380:    GOSUB  213B0:    NORMAL    :    PR I NT   :    PRINT   :  PRINT 

Hfo  SA  I  SAA/*K,?^399):L°  =   10:HI   =   150001   G0SUB  39300:  SA  -  CV 
2040  A  =   11.6:B  =   1. 18 
2050  C  =  3. 86: E  =  2.32 
2060  EX   =  .94 
2070  RF  =  100 


2080  X   =  SA  *   XI    *    (RF  -    (AD)    *   31    *   24   *   3600   /    <L1    *    10  -  9) 
2090  X  =  X   *    (A  +    (B  *  TH)    +    (C  *  Tl)    -   (E  *  AD)    /    (RF  -  (At! 
2100  Y  =  Yl   *  EX  *   II   *   10  *  6  *  31   *  SA  /    (LI   *   10  ~  9) 
2110  PI   =    (i   029  *  Y)    -    (.065  *   X)    -    (.245  *  Y  -  2)    +    (  0018  *   X  '  2) 
+    (.0215  *   Y       3)  Z) 


2120  X  =  SA  *   XI    *    (RF  -    (A2))    *   28  *   24   *   3600   /    (L2  *   10  -  9) 

2130  X  =  X   *    (A  ♦    (B  *   TH)    +    (C   *   T2)    -    <E   *  A2>)    /    (RF  - 

2140  Y  -  Yl   *  EX   *   12  *   10  -  6  *  28  *  SA  /    (L5J   10  '9)  ' 

P2  7.02i5,2I  *Y  *   X>    "    ("245  *  V  -  2)   +  (.0018 

2160  X   =  SA  *   XI    *    (RF  -    (A3))    *   31    *  24   *   3600   /    (L^  *   10  - 

2 III  y :  JiVJxYkY?; :  ic  * T3)  - (E  * A3>> ' <R*  -  <a3> 

^1BU  Y  m  Yl   *  EX  *  13  *   10    s  6  *  31   *  SA  /    <L*  *   10  '-  9) 

2190  pl  T.wiS2!  *      ('0£>5  * X)  ~ (-245  * Y  * 2)  +  (-0018  * x  *  2> 

I?00  5  =  ®A  *  *'   *  -    (A4))    *  30  *  24  *  3600  /    (L4  *   10  a  9) 

2210  X  =  X  *    (A  ♦   (B  *  TH)   +   (C  *  T4)   -   (E  *  A4)  )   /    (RF  -  (A4) 
222U  Y  =  Yl   *  EX   *   14  *   10  -  6  *  30  *  SA  /    (L4  *   10  -  9) 

!  7.02l52I  *  "z~>  (-°65  *  X>  "  ('245  *  v  ^  2)  +  (-°°18  *  X  *  2> 

™°  J  =  *  "    <A5>>    *  31    *  24   *   3600  /    (L5   *   10  -  9) 

"J|J  x  =  X  *  (A  +  (B  *  TH)  +  (C  *  T5)  -  (E  *  A5)  )  /  (RF  -  /flc/, 
-260  Y  -  Yl    *   EX   *    15  *    10  ~  6   *   31    *   SA  /    (L5  *    10    '   9)        ^ ' 

+   T.02iS2I  J  -}3;    ('°65   *   X>    "    (-245  *   Y  ""   2>    +    (-°018  *   X  -  2) 

2280  X  =  SA  *   XI    *    (RF  -    (A6) )    *   30  *   24   *  3600   /    (LI    *    10  -  9) 

229<,  X   =   X   *    (A  +    (B  *   TH)    ♦    (C   *   T6)    -    (E  *  A6)  )    /    (RF  -  <A6>) 

2300  Y  -  Yl   *  EX  *   16  *   10  -  6  *  30  *  SA  /  (L6  *   10  -  9) 

t  ;.o2i52l  *  i}3r.      * X) " (-245  * Y  ~  2>  +  ('0018  *  X  ^  2> 

2320  X   =  SA   *    XI    *    (RF  -    (A7) )    *   31    *  24   *   3600  /    (L7   *    10  *  9) 
2330  X  =  X   *    (A  ♦    (B  *  TH)    +    (C  *   T7)    -   (E   *  A7)  )    /    (RF  -  (A7> 
2340  Y  -  Yl  *  EX  *  17  *  10  -  6  *  31  *  SA  /    (L7  *   10  -9) 

I  T.02l52I  V  (-°65  *   X>    "   ("245  *  Y       2>    +    (-0018  *  X  -  2) 

2360  X  =  SA   *   XI    *    (RF  -    (A8) )    *   31    *   24   *   ^600  /    (L8  *    10  ~  9) 

™'       =   X   *    <A  +    <B   *   TH)    +    (C   *   T8)    -    <E   *  AS))    /    (RF  -  (AS)) 

2.380  Y  =  Yl    *   EX   *    18  *    10  -  6   *   31    *   SA   /    (L8  *    10  -  9) 

2390  P8  =    (1.029   *   Y)    -    (.065   *    X)    -    (.245   *   Y  --  2)    +    (.0018  *   X   -  2) 

+    (.0215   *   Y  -•■  3) 
2400  X   =  SA  *   XI    *    (RF  -    (A9))    *   30   *   24   *   3600   /    (L9  *    10  ■'-  9) 
2410  X  =   X   *    <A  +    (B  *   TH)    +    (C   *   T9)    -    (E   *  A9) )    /    (RF  -    (A9) ) 
2420  Y  =  Yl   *  EX  *   19  *   10  -  6  *  30  *  SA  /    (L9  *   10  -  9) 
2430  P9  =    (1.029  *  Y>    -    (.065  *   X)    -   (.245  *  Y  -  2)    +   (.0018  *   X  '-  2) 

+    (.0215  *  Y  ••■  3) 
2440  X  =  SA   *   XI    *    (RF  -    (AO))    *   31    *   24   *   3600   /    (LO  *    10  ~  9) 
2450  X  =   X   *    (A  +    (B  *   TH)    +    <C   *   TO)    -    <E   *  AO))    /    (RF  -  (AO)) 
2460  Y  =  Yl    *  EX   *    10  *    10  •    6   *   31    *   SA.  /    (LO  *   10  -  9) 


2470  P0  =    (1.029   *   Y)    -    (.065   *    X)    -    (.245   *   Y       2)    +    (.0018   *   X  2) 

+    (.0215   *   Y  •-•  3) 
2480   X   -  SA   *   XI    *    (RF  -    (AY))    *   30   *   24   *   3600   /    (LA   *    10  ••  9) 
2490   X   =   X    *    (A  +    (B   *   TH)    +    (C   *   TA)    -    (E   *   AY))    /    (RF  -  (AY)) 
2500  Y  =  Yl    *   EX   *    IA   *    10    -  6   *   30   *   SA   /    (LA   *    10  "  9) 
2510  PA  =    (1.029   *   Y)    -    (.065   *   X)    -    (.245   *   Y  rt  2)    +    (.0018   *    X   A  2) 

+    (.0215   *   Y  •-•  3) 
2520  X   =  SA  *    XI    *    (RF  -    (AZ))    *   31    *  24   *   3600   /    (LB   *    10  *  9) 
2530  X   =   X   *    (A  +    (B   *   TH)    +    (C   *   TB)    -    (E   *   AZ))    /    (RF  -  (AZ)) 
2540  Y  =  Yl    *  EX   *    IB   *    10  •"    6   *   31    *  SA   /    (LB   *    10  "  9) 
2550  PB  =    (1.029   *   Y)    -    (.065   *    X)    -    (.245   *   Y  •"•  2)    +    (.0018  *    X   '"•  2) 

+    (.0215   *   Y  ■-  3) 
2560      IF  PI    >   1   THEN  PI   =  1 
2570     IF  FA  =   1    THEN  PI   =  PI    *  .8 
2580      IF   FA  -   3   THEN  PI    =   PI    *  .9 
2590      IF  P2   >   1    THEN  P2  =  1 
2600      IF  FA  =   1    THEN  P2  =  P2  *  .8 
2610      IF  FA  =  3  THEN  P2  w  P2  *  .9 
2620      IF   P3    >    1    THEN   P3   =  1 
2630      IF   FA  =    1    THEN  P3   =   P3   *  .8 
2640      IF  FA  =  3  THEN  P3  =  P3   *  .9 
2650      IF  P4   >   1   THEN  P4  =  1 
2660     IF  FA  =   1   THEN  P4   =  P4   *  .8 
2670     IF  FA  =  3  THEN  P4  «  P4  *  .9 
2680      IF  P5   >   1   THEN  P5  =  1 
2690     IF  FA  =   1   THEN  P5  »  P5  *   . 8 
2700     IF  FA  =  3  THEN  P5  -PS.*  -9 
2710     IF  P6   >   1   THEN  P6  ■  1 
2720     IF  FA  =   1   THEN  P6  =  P6  *  .8 
2730     IF  FA  =  3  THEN  P6  =  P6  *  .9 
2740     IF  P7   >   1   THEN  P7  =  1 
2750     IF  FA  =   1   THEN  P7  -  P7  *  .8 
2760     IF  FA  =  3  THEN  P7  =  P7   *  .9 
2770     IF  P8   >   1   THEN  P8  =  1 
2780     IF  FA  =    1   THEN  P8  =  P8   *  .8 
2790      IF  FA  =  3  THEN  P8  =  P8   *    .  9 
2800      IF  P9   >   1   THEN  P9  =  1 
2810      IF  FA  =   1   THEN  P9  =  P9   *  .8 
2820     IF  FA  =  3  THEN  P9  =  P9  *    „ 9 
2830      IF   PO   >    1    THEN  PO   =  1 
2840      IF  FA  -    1   THEN  PO  =  PO   *  .8 
2850      IF  FA  =  3  THEN  PO  =  PO   *  .9 
2860      IF  PA   >    1   THEN  PA  =  1 
2870     IF  FA  =   1   THEN  PA  =  PA   *    . 8 
2880      IF   FA   =   3   THEN  PA   =   PA   *  .9 
2890      IF   PB    >    1    THEN  PB  =  1 
2900     IF  FA  =   1   THEN  PB  =  PB   *  .8 
2910      IF   FA  =   3   THEN  PB   =   PB   *  .9 

2920  LQ  =   LI    +  L2   +   L3   +   L4   +   L5   +   L6   +   L7   +   L8   +   L9  +  LO   +  LA   +  LB 
2930  QQ  =    (PI    *  LI)    +    <P2   *   L2)    +    (P3   *   L3)    +    <P4   *   L4)    +    (P5   *   L5)  + 
<P6   *   L6)    +    <P7   *   L7)    +    (P8   *   L8)    +    (P9   *   L9)    +    (PO   *   LO)    +    (PA  * 
LA)    +    (PB   *  LB) 
2940  AS  =  QD   /  LQ 
2950      IF   AS    >    1    THEN  AS   =  1 

2960     HOME   :    INVERSE   :    GOSUB  21380:    PRINT   "  *   FRACTION  OF  HOT  WATE 

R";:    GOSUB  21380:    PRINT    "  SUPPLIED  BY  SOLAR"::    GOSUB  21380 

:    GOSUB  21380:  NORMAL 

2965  GOSUB  22200:  TM  -  Y5  +   1:    GOSUB  22100:  PRINT 

2966  LM  =  5:    GOSUB  22120:  PRINT 

2970     PRINT  JA$;"    =";    TAB (    20);    INT    (PI    *    100)    /  100 


2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3085 
3090 
3 1 00 
3110 
3115 

3120 
3130 
3140 
3150 
3160 
3165 
3 1 80 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 


3280 

3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 


3500 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 
LM  =  0 

PRINT 

LET  NS  = 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 
ZR  - 


FB* ;  " 

= 11  • 

TAB  ( 

20)  ; 

INT 

<P2 

* 

100) 

/ 

1 00 

MR*;  " 

=  "  « 

TAB  < 

20)  ; 

INT 

<P3 

* 

100) 

/ 

100 

AP*;  " 

=  »  • 

TAB  ( 

20)  ; 

INT 

<P4 

* 

100) 

/ 

100 

MY* ;  " 

=  "  ■ 

TAB  ( 

20)  ; 

INT 

(PS 

* 

1 00 ) 

/ 

100 

JE*;  " 

=  "  . 

TAB  ( 

20)  ; 

INT 

<P6 

* 

100) 

/ 

100 

JL*;  " 

—  "  . 

TAB  < 

20 )  ; 

I  NT 

(P7 

* 

100) 

/ 

1 00 

AG* ;  " 

=  "  . 

p 

TAB  ( 

20)  ; 

INT 

<P8 

* 

100) 

/ 

100 

SP*;  " 

_  II  . 

TAB  ( 

20)  ; 

INT 

<P9 

* 

1 00 ) 

/ 

100 

0T*;  " 

—  "  . 

TAB  < 

20)  ; 

INT 

<P0 

100) 

/ 

1 00 

NV*;  " 

_  II  . 

TAB  ( 

20)  ; 

INT 

(PA 

* 

100) 

/ 

1 00 

DC*;  " 

—  II  . 

f 

TAB  ( 

20)  ; 

INT 

<PB 

* 

100) 

/ 

1 00 

GOSUB  22120 
PRINT    "ANNUAL  ="; 
(1   -  AS)    *  (AG). 


TAB(   20);    INT    (AS   *    100)    /  100 


"HIT  ANY  KEY  TO  CONTINUE 


GET  AK* :    HOME   :  PRINT 


PRINT 
PRINT 
TEXT  : 
IF  FA 


"ANNUAL  KWH  OF  ELEC  USED' 
"AS  BACKUP  FOR  SOLAR  =  "; 
INT    (ZL   *  AS   *  10) 


TAB (   20);    INT  (NS) 


:  PRINT 
:  PRINT 
HOME 


"SOLAR  HEAT=  "ZR  /  10"  MBTUS " 
"PRESS  ANY  KEY  TO  CONTINUE  . 


GET  AK* 


FA 
FA 
FA 
FA 
FA 


1 

THEN 

FC 

500 

1 

THEN 

VC 

215 

THEN 

FC 

1200 

THEN 

VC 

325 

3 

THEN 

FC 

700 

3 

THEN 

VC 

325 

-  TX) 

/ 

1 00 

FC  + 

(VC 

* 

SA) 

IF 
IF 
IF 
IF 
IF 

KX  =  (100 
LET  CS  = 
GOTO  3300 

HOME   :    INVERSE   :    GOSUB  21380: 
T  OF";:    GOSUB  21380:    PRINT  " 
80:    GOSUB  21380:    NORMAL   :  PRINT 

PRINT   "ENTER  THE   INTIAL  COST  OF 
R-LEAVE  OFF  *  SIGNS  AND  COMMAS. " 
CV  =  CS:LO  =  0| HI    =  50000:    GOSUB  39300: CS 
OA        (CS  -   1000)    *  .025 
OB  =   1000   *  .05 
SX   =  OA  +  QB 
IF  CS  <    1001    THEN  SX   =  CS   *  .05 


PRINT    "        *   ESTIMATED   INITIAL  COS 
INSTALLED  SYSTEM.*";:    GOSUB  21^ 

PRINT   :  PRINT 
THE  SYSTEM  TO  THE  NEAREST  DOLLA 


CV 


IF 


CI   =  (CS 
OM  =  CS  * 
SI   =  CG  * 
N  =  20 

IF   ZQ  -  1 

IF  ZQ 
QB 
RR 
DP 
ZZ 
DP 


125  THEN  SX  =  125 
*   KX)    -  SX 
.  01 

(AG  -  NS) 


THEN  D  = 
=  2  THEN  D  = 

=  D 

=    (RQ  +  ES)    /  2 
*»   DP   +  1 

=      LOG    ( ( 1    +  RR) 
=     LOG    (1    +  CI  / 
IF  AA  =  2  THEN  D  m 
PRINT 

HOME   :    INVERSE   :    GOSUB  21380:    PRINT    "  *   EXPECTED  FIRST  YEAR 

";:    GOSUB  21380:    PRINT    "  ECONOMIC  PERFORMANCE"; :    GOSUB  21 

380:    GOSUB  21380:    NORMAL    :    PRINT    :    PRINT   :  PRINT 

PRINT   "CAPITAL  CST   OF   SYSTM  =";    TAB (   25);"*";    INT  (CS) 


.  046 
.  0278 


/    (1    +  D)  ) 
<S1    -  OM)  * 
.  0278 


(RR 


-  D) )    /  ZZ 


3510  RX  =   1   -  KX 
3520  FX   =  CS   *  RX 
3530  PRINT 
3540  PRINT 

3550  PRINT 

3560  PRINT 

3570  PRINT 

3580  PRINT 

3590  PRINT 

3600  PRINT 

3610  HOME 

:    PRINT  " 

:  NORMAL  :  PRINT  : 
PRINT  "DISCOUNTED 
PRINT 


"FEDERAL  TAX   CR.=";    TAB (    25);"*";  INT 
•    PRINT    " STATE  TAX   CR.=":    TAB (  25);"*" 
"FIRST   YR  CST  W/TX   CR.=";    TAB < 
"0  8<  M  COSTS  TAB  (  25);"*"; 

"FIRST   YR  SAVINGS  =";    TAB  <  25) 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT  ' 
INVERSE 


(FX) 

;    INT  (SX) 
25);"*";  INT 
INT  (OM) 

;"*";    INT  (SI) 


(CI 


KEY 


3620 

3630 

364  0 

3650 

3660 

36  70 

3680 

3690 

3700 

37 1 0 

3720 

3730 

3740 
3760 
3770 
3780 
3790 
3800 
3810 
3820 
3S30 
3850 
3860 
3870 
3880 
3890 
3900 
3910 
3920 
3930 
3940 
3950 
3960 
3970 
3980 

3990 
4000 
4010 
4020 

4030 
4040 
4050 
4060 

4070 


HIT  ANY 
:  PRINT 
ECONOMIC 
PR  I  NT 
PAYBACK 


TO  CONTINUE   .    .    .";:    GET  AK* 
*   EXPECTED  LONG  TERM";:    GOSUB  21 
PERFORMANCE";:    GOSUB  21380:    GOSUB  21 


INT    (DP)"  YRS" 


S  =  O 
RP  =  0 
PV  =    ( ( 1 
M  =  OM  * 


RG 
NR 
NM 
EP 
EN 
NC 
LC 


O 
0 

200 
(NS 
(AG 
NR  ■ 


+  D) 
PV 


*  .01 

*  CG) 

*  CG) 

*  S  + 


N 


1)    /    (D   *    (1    +  D)  N) 


=  CI   -  S  + 


* 

* 

NM 
M 


PV 
PW 
PW 

+  RG 
h  RP  - 


'LI FECYCLE  COST 


h  EN 
EP 

/SOLAR 


INT    (NC) :NC  =  NC   /  100 


PRINT 
PRINT 
NC  =  NC   *    100: NC 
PRINT 

PRINT   "LI FECYCLE   COST  /ELECTRICITY 
PRINT 
NC  =  NC  -  LC 
PRINT 

PRINT   " NET  SAVINGS   

GOTO  5150 
SA  =  SA  +  1 
ZB  =  ZB  +  1 

NORMAL 
XX   =  NS 

A  =  11.6:B  =  1.  18 
C  =  3.86:E  =  2.32 
EX  =  .94 

RF  =  100 

X  =  SA  *  XI  *  (RF  -  (Al) )  *  31  * 
X   =   X*(A+(B*TH)    +  (C*T1) 


*";    INT  (LC) 


$";    INT  (NC) 


*";    INT  (NC) 


Y  =  Yl  *  EX  * 
PI    =    (1.029  * 

+  (.0215  *  Y 
X  —  SA  *  X 1  * 
X   =   X    *    (A  + 

Y  =  Yl  *  EX  * 
P2  =    ( 1 . 029  * 

+  (.0215  *  Y 
X  =  SA  *  XI  * 
X   =   X    *    (A  + 

Y  =  Yl  *  EX  * 
P3   =    (1.029  * 

+  ( .  02 1 5  *  Y 
X   =   SA   *    XI  * 


10  "•  6  *  31  * 
(.065   *    X)  - 


I  1  * 

Y)  - 

"•  3) 

(RF  -  (A2) )  *  28  * 
(B   *    TH)    +    (C    *  T2) 

12  *  10  6  *  28  * 
Y )    -    ( . 065   *    X )  - 

"  3) 

(RF  -  ( A3) )  *  31  * 
(B    *   TH)    +    (C   *  T3) 

13  *  10  ''-  6  *  31  * 
Y )    -    < . 065   *   X  >  - 

•-■  3) 

(RF   -    ( A4) )    *   30  * 


24    *   3600   /    (LI    *  10 
-   (E  *  Al )  )    /    mr  - 
SA   /    (LI    *    10  -  9) 
(.245   *   Y  •-    2)    +  (.0018 


■  9) 
(Al)  ) 


*  X 


24    *  3600   /    (L2   *    JO    -  9) 

-  (E  *   A2)  )    /    (RF   -    (A2)  ) 
SA   /  (L2   *    10    -  9) 

(.245  *    Y   ■-■   2)    +    (.0018   *  X 

24    *  3600   /    (L3    *    10   ■""  9) 

-  (E  *   A3)  )    /    (RF   -    (A3)  ) 
SA   /  (L3   *    10  -  9) 

(.245   *   Y     ■  2)    +    (.0018   *  X 


24   *   3600   /    (L4    *  10 


9) 


4080 
4096 
'4100 

4110 
4 1  20 
4130 
4140 

4 1 50 
4  1 60 
4 1 70 
41  BO 

4190 
4200 
4210 
4220 

4230 
4240 
4250 
4260 

4270 
4280 
4290 
4300 

4310 
4320 
4330 
4340 

4350 
4360 
4370 
4380 

4390 
4400 
4410 
4420 

4430 
4440 
4450 
4460 
4470 
4  480 
44  90 
4  500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 


X  =  X   *    (A  +    (B  *  TH)    +    (C  *  T4)    -  (E 

Y  =  Yl  *  EX  *  14  *  10  '  6  *  30  *  SA  / 
P4  =    (1.029   *   Y)    -    (.065   *    X)    -  (.245 

+    (.0215  *   Y  3) 
X  =  BA  *   XI   *    (RF  -    (A5>)    *  31    *  24  * 
X  =  X   *    (A  +    (B  *   TH)    +    (C   *   T5)    -  (E 

Y  =  Yl  *  EX  *  15  *  10  •-  6  *  31  *  SA  / 
P5  =    (1.029  *   Y)    -    (.065  *   X)    -  (.245 

+    < .  02 15  *   Y       3 ) 
X   =  SA   *   XI    *    (RF   -    <A6>)    *   30   *   24  * 
X  =  X   *    (A  +    (B  *   TH)    +    (C  *  T6)    -  (E 

Y  =  Yl  *  EX  *  16  *  10  6  *  30  *  SA  / 
P6  =    (1.029   *   Y)    -    (.065   *    X)    -  (.245 

+    (.0215  *  Y  -"■  3) 
X  m SA  *  XI   *    (RF  -    (A7)>    *  31   *  24  * 
X   =  X   *    (A  +    (B  *   TH)    +    (C   *   T7)    -  (E 

Y  =  Yl  *  EX  *  17  *  10  ••  6  *  31  *  SA  / 
P7  =    (1.029   *   Y)    -    (.065   *    X)    -  (.245 

+    (.0215  *  Y  3) 
X   =  SA   *   XI    *    (RF   -    (AB))    *   31    *   24  * 
X   =   X   *    (A  +    (B   *   TH)    +    (C   *   T8)    -  (E 

Y  =  Yl  *  EX  *  18  *  10  '•  6  *  31  *  SA  / 
PB  =    (1.029   *   Y)    -    (.065   *    X)    -  (.245 

+    ( . 02 1 5  *   Y  '    3 ) 
X   =  SA   *   XI    *    (RF  -    (A9>)    *   30   *  24  * 
X   =  X   *    (A  +    (B   *   TH)    +    (C   *   T9)    -  (E 

Y  =  Yl  *  EX  *  19  *  10  ~  6  *  30  *  SA  / 
P9  =    (1.029   *   Y)    -    (.065   *    X)    -  (.245 

+    (.0215  *  Y  '  3) 
X  =  SA  *  XI   *    (RF  -    (AO) )    *  31   *  24  * 
X  =  X   *    (A  +    (B  *  TH)    +    (C   *  TO)    -  (E 

Y  -  Yl  *  EX  *  io  *  10  6  *  31  *  SA  / 
PO  =    (1.029   *   Y)    -    (.065   *   X)    -  (.245 

+    (.0215  *  Y  -  3) 
X   =  SA  *   XI    *    (RF  -    (AY))    *   30   *   24  * 
X  =  X   *    (A  +    (B  *  TH)    +    (C   *  TA)    -  (E 

Y  =  Yl  *  EX  *  IA  *  10  •-•  6  *  30  *  SA  / 
PA  =    (1.029   *   Y)    -    (.065  *    X)    -  (.245 

+    (.0215  *   Y  '"•  3) 
X   =  SA   *   XI    *    (RF  -    (AZ))    *   31    *   24  * 
X   =  X    *    (A  +    (B   *   TH)    +    (C   *   TB)    -  <E 

Y  =  Yl  *  EX  *  IB  *  10  A-  6  *  31  *  SA  / 
PB  =    (1.029   *   Y)    -    (.065  *    X)    -  (.245 

+    (.0215   *   Y  3) 


*  A4) )    /    (RF  -    (A4) ) 
(L4  *   10  '-  9) 

*  Y  -  2)    +    (.0018  *   X  ^  2) 

3600   /    (L5   *    10  •-  9) 

*  A5) )    /    (RF  -    (A5) ) 
<L5  *    10  ••••  9) 

*  Y  ■--  2)    +    (.0018  *   X  ■'-  2) 

3600  /    (LI   *   10  •'  9) 

*  A6> )    /    (RF  -    (A6) ) 
(L6   *    10  9) 

*  Y  -"■  2)    +    (.0018  *   X   •'-  2) 

3600  /    <L7  *   10  •■-  9) 

*  A7) )    /    (RF  -    (A7)  ) 
(L7   *    10  A  9) 

*  Y       2)   +    (.0018  *   X  •'-  2) 

3600  /    (L8   *    10  9) 

*  AB) )    /    (RF  -    (A8) ) 
(L8  *    10  '"-  9) 

*  Y  -■-  2)    +    (.0018   *   X  2) 

3600  /    (L9   *   10  *  9) 

*  A9) )    /    (RF  -    (A9)  ) 
(L9  *    10  •'-  9) 

*  Y  '-■  2)    +    (.0018  *  X  *  2) 

3600  /    (LO  *   10  A  9) 

*  AO) )    /    (RF  -    (AO) ) 
(LO  *   10  *  9) 

*  Y  *  2)   +   (.0018  *   X  ^  2) 

3600  /    (LA  *   10  *  9) 

*  AY) )    /    (RF  -    (AY) ) 
(LA  *    10  •■•  9) 

*  Y  '-•  2)    +    (.0018  *  X  "  2) 

3600   /    (LB   *    10  9) 

*  AZ) )    /    (RF  -    (AZ) ) 
(LB  *   10  9) 

*  Y       2)   +    (.0018  *   X  2) 


IF  PI  >  1  THEN  PI 
IF  FA  =  1  THEN  PI 
IF  P2  >  1 
1 
1 
1 
1 
1 
1 


THEN  P2 
THEN  P2 
THEN  P3 
THEN  P3 
THEN  P4 
THEN  F4 
THEN  P5 
1   THEN  PS 
1   THEN  P6 
1    THEN  P6 
1   THEN  P7 
IF  FA  =   1   THEN  P7 
IF  P8   >   1    THEN  P8 


IF  FA  = 
IF  P3  > 
IF  FA  = 
IF  P4  > 
IF  FA  f 
IF  P5  > 
IF  FA  = 
IF  P6  > 
IF  FA  = 
IF  P7  > 


1 

PI 
1 

P2 
1 

P3 
1 

P4 
1 

PS 
1 

P6 
1 

P7 
1 


*  .8:  IF  FA  = 

*  .8:  IF  FA  = 

*  .8:  IF  FA  = 

*  .8:  IF  FA  =  : 

*  .8:  IF  FA  = 

*  .8:  IF  FA  = 

*  .8:  IF  FA  = 


IF  FA  =   1   THEN  PB  =  P8   *    .8:    IF  FA 


3  THEN  PI 
3  THEN  P2 
3  THEN  P3 
3  THEN  P4 
3  THEN  PS 
3  THEN  P6 
3  THEN  P7 
3  THEN  PB 


PI    *  .9 

P2  *  .9 

P3   *  .9 

P4    *  .9 

PS   *  .9 

P6   *  .9 

P7   *  .9 

PB   *  .9 


4590 

I  F 

F'9  > 

1 

THEN 

P9 

= 

1 

P9 

4600 

IF 

FA  = 

1 

THEN 

P9 

= 

P9 

* 

.  8: 

IF 

FA  = 

3 

THEN 

P9  - 

* 

.9 

461  0 

IF 

P0  > 

1 

THEN 

PO 

= 

1 

4620 

IF 

FA  = 

1 

THEN 

PO 

PO 

.8: 

IF 

FA  - 

3 

THEN 

PO  = 

PO 

* 

.  9 

4630 

IF 

FA  > 

1 

THEN 

PA 

= 

1 

4640 

IF 

FA  = 

1 

THEN 

PA 

PA 

.  8: 

IF 

FA  = 

3 

THEN 

PA  = 

PA 

* 

.9 

4650 

IF 

PB  > 

1 

THEN 

PB 

1 

4660 

IF 

FA  = 

1 

THEN 

PB 

PB 

* 

.8: 

IF 

FA  ■ 

3 

THEN 

PB  = 

PB 

* 

.  9 

4670  LO  =  LI    +  L2  +  L3  +  L4  +  L5  +  L6  +  L7  +  L8  +  L9  +  LO  +  LA  +  LB 
46BO  00  =    (PI    *  LI)    +    (F2   *   L2)    +    <P3   *   L3)    +    <P4   *   L4)    +    (P5  *   L5)  + 
<P6   *   L6)    +    <P7   *   L7)    +    (PB   *   L8)    +    <P9   *   L9)    +    (PO   *   LO)    +    (PA  * 
LA)    +    (PB   *  LB) 
4690  AS  =   00   /  LO 
4700      IF   AS    >    1    THEN  AS   -  1 
4710   SS   =    (1    -    AS)    *  AG 
4740      IF   FA  =    1    THEN  FC   =  500 
4750      IF  FA  =    1    THEN  VC  =  215 
4760      IF   FA  =   2   THEN  FC   =  1200 
4  770      IF    FA   =   2   THEN  VC   =  325 
4780      IF   FA  =   3   THEN  FC   =  700 
4790      IF  FA  =  3  THEN  VC  =  325 
4800  CS  =  FC  +    (VC   *  SA) 
4810  KX   =    (100  -   TX)    /  100 
4820  OA  =    (CS  -    1000)    *  .025 
4830  OB  =    1000   *  .05 
4840  SX   =  OA  +  OB 

4850      IF  CS  <    1000  THEN  SX   =  CS   *  .05 
4860      IF   SX    >   125  THEN  SX   =  125 
4870  CI   =    (CS   *   KX)    -  SX 
4880  ON  =  CS   *  .01 
4890  N  =  20 

4900      IF   ZQ  =   1   THEN  D  =  .046 
4910      IF   ZO  =  2  THEN  D  =  .0278 

4920  S  =  0 
4930   RP   =  0 

4940  PV  =    ((1   +  B)    -    N-1)/(D*(1+B)    A  N) 
4950  M  =  DM   *  PV 
4960  RG  =  0 
4970  NR  =  0 

4980  NM  =   200   *    .01    *  PV 

4990   EP   =    (SS   *   CG)    *  PW 

5000   EN  =    (AG   *   CG)    *  PW 

5010  NC  =  NR  -  S  +  Nti  +  RG  +  EN 

5020  LC  =  CI   -  S  +  M  +  RP  +  EP 

5030  NS   =   NC   -  LC 

5040  NORMAL 

5050  PRINT 

5060     PRINT    "  " 

5070     PRINT    " SQR  FT-   "    INT    (SA   *  10.76) 

5080     PRINT    " SOLAR  FRACTION-    ";(    I NT    (AS   *    100))    /  100 
5090     PRINT    "NET   SAVINGS^  *";    INT  (NS) 

5100     PRINT    "  " 

5110      IF   NS    >   XX   GOTO  3870 
5120      IF   NS   =   0   THEN  3870 
5130      IF   ZB   =    1    THEN  3870 
5140  RETURN 

5150  SIR   =    ((EN   -   EP)    -    ( M   -   NM )  )    /    CI    -    CS  +   R > 

5160  PRINT 

5170  SN  =   SIR   *  100 

5180  SN   =      I  NT  (SN) 


519.0 
,•5200 
5210 


5240 
5250 
5260 
5270 


5290 
5300 
5310 
5320 
5330 
5340 
5350 
5360 
5370 
5380 
5390 

5400 
54 1 0 

5420 
5430 
5435 
5440 

5450 
5460 
5470 
5480 


5490 
5500 
5510 
5520 
5530 
5540 


5550 
5560 
5570 
5580 
5590 
5600 
5610 
5620 
5630 
5640 
5650 
5660 


SN  =  SN  /  100 

PRINT   :    PRINT    "SAVINGS   INVESTMENT  RAT T n  .-  r-h, 

PRINT   :   PRINT   :    PRINT  "    "  * 

PRINT    "HIT  ANY  KEY  TO  CONTINUE   .    .  GEJ 

eES?  loilS^,^2^'  ™"  "  '       1  SSld  you  l|ke  to  s 

21380:    PRINT   "  USED   IN  THE  Inai        ^IAL  ASSUMPTIONS";  i  GOSUB 

21380:    NORMAL   :    PRINT   "  PRINT     "pRlT  ''  2138°!  B°SUB 

PRINT   "1)    YES":    PRINT   "2)    NO":  PRINT 

HOME^:L°  =   1:HI   =  25    G0SUB  393001    IF  CV  =  2  THEN  5400 

^■i^ '»sk  ss'SE  nor^c;app~ters  for  « 

PRINT 

PRINT  TAB<  12);"DI8C0UNT  RATE=     "; (D  +  .08)    *  100"%" 

PRINT  TAB (  12);"C0ST   OF  ELEC. =  *"CG"/KWH" 

PRINT"  TAE(  12>i  "FUEL  PRICE  ESCAL  RATE=   "<ES  *    100)    *  B»X- 

PRINT     TAB  (    12)  ;  "INFLATION  RATE  =  87." 
PRINT 

PRINT     TAB (    12);  "LIFE  OF   INVESTMENT"     20  YFAPft» 

PRINT   :    PRINT   ,    PRINT    "HIT  ANY  n  'coNTINul  GET  AK$ 

HOME  :    IF  PTt  =  "N«  THEN  5440 

PRINT     1)    YES":   PRINT   "2)    NO":  PRINT 
pV0tvL°  c>   ltHI   =  2!    G0SUB  393005    IF  CV  ■   1    THEN     GOSUB  5960 
^::^^°™^r^T   "      *  WHERE  TO  NOW?" ;  :  GOSUB 
Cphn  °:L0  =  °:HI   =   128:    GOSUB  39300 

CV3f  XCi?8^  21380:    N°RMAL   =    ^  '   ^  I  NT  r  PR  I  NT      <FR-UL>?"'=  G0S^ 
NORMAL  =   -°5:HI    =  GOSUB  39300:  X(l)    =  CV 

IF  XU)    -  0  THEN  X(l)    =  .68 
Yd)   =  0 

V  i-V^i^lH^fS^T^,;^^  PRINT   ,  PRINT™ 
XI   =  0:  Yl    =  0 
XY  =  0:X2  =  0 

FOR   1=1   TO  N 
XI   -   XI   4-   X  ( I ) 
VI    -  Yl   +  Y(I) 
XY  =  XY  +  XU)    *  Y(I) 
X2  =  X2  +  X(I)    *  X<I) 

NEXT  I 
J   =  N   *   X2  -   XI    *  XI 

IF  J   <      >  0  THEN  5670 

PRINT   "NO  SOLUTION  FOUND":    FOR  X   =   1    TO   1000:    NEXT   X : X  ( 1 )    =  0:Y 


5670 
5680 
5690 
5700 
57 1 0 
5720 
5750 
5760 
5770 
5760 
5790 
5800 
5810 
5820 
5830 
5840 
5850 
5860 
5870 
5880 
5890 
5900 
5910 
5920 
5930 
5940 
5950 
5960 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 


(2) 
M  = 
M  « 
B  = 
IF 

RETURN 
REM 
K  =  ES 
P  =  0 
E   =  .06 
RQ  =  E 
AD  = 
AB  = 
FDR 
P  (N> 


=  0:  GOTO 
<N  *   XY  - 
INT  (1E3 
(Yl    *  X2 
ABS  (B) 


55 1 0 

XI    *   Yl)    /  J 
*   M  +   .5)    /  1E3 
•    XI    *   XY)    /  J 


=  B  THEN  B* 


m  +  .i 


D) 
D) 
TO 
(  ( 1 


/ 
/ 

10 


<  1 
(  1 


AO) 


11   TD  20 

/    ( ( 1   +  AB ) 


1  TO 

P(N) 


10 


20 


<  ( 1 

<  ( 1 
N  = 
=  1 

NEXT  N 
FOR  N  « 
V(N)    =  1 
NEXT  N 
FOR  N  = 
P  =  P  + 
NEXT  N 

FOR  N  =   11  TO 
V  =  V  +  V(N) 

NEXT  N 
PW  =  P  +  V 
PZ  =  PW 

RETURN 

REM     PRINT  SUBROUTINE 
D*  =   " " 
PRINT  D$ ; "PR#1 " 
PR#  1 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


E)  ) 
K )  ) 


N) 


N) 


****** 
PRINT 
PR  I  NT 
PRINT 


SOLAR   HOT   WATER   HEATING   ANALYSIS  ******" 


'LOCATION=    "RE*",  MONTANA" 
"NO.    PERSONS   IN  HOUSEHOLD"  ";ZP 


PRINT    "HOT   WATER   USED   PER  DAY 


INT    (Q   *    ZP) "  GALLONS 


6050  PRINT 

6060  PRINT   " HOT  WATER  TEMP. ■   ";<TH   *    1.8)    +  32"   DEGREES " 

6070  PRINT 

6080  PRINT   "ANN.    KWI!  USED  FOR  HOT  WATER"   " ;    INT  (AG) 

6090  PRINT 

6100  PRINT   "OPTIMAL  FT2  OF  COLLECTOR"  INT    (SB   *  10.7639) 

6110  PRINT 

6120  PRINT   "AREA  ANALYZED" 


6130  PRINT 
6140  PRINT 


INT    (SA   *    10. 764) "   FT2 " 
ANN.    SOLAR   FRACTION   OF   HOT   WATER-    "    INT    (AS   *    100)    /  100 


6150  PRINT 

6160  PRINT    "ANN.    SOLAR  ENERGY  FROM  SYSTEM"   " ; ZR   /    10"   MBTUS " 

6170  PRINT 

6180  PRINT   "ANN.    KWH  USED  FOR  BACKUP"   " ;    INT    (NS   *    100)    /  100 

6190  PRINT 

6200  PRINT   "INTIAL  COST  OF  SYSTEM"  *";    INT  (CS) 

6210  PRINT 

6220  PRINT   "FEDERAL  TAX   CREDI T=  INT  (FX) 

6230  PRINT 

6240  PRINT   "STATE  TAX   CREDIT"  INT  (SX) 

6250  PRINT 


6260 

PR  I  NT 

6270 

PRINT 

6280 

PRINT 

6290 

PRINT 

6300 

PRINT 

6310 

PRINT 

6320 

PRINT 

6330 

PRINT 

6340 

PRINT 

6350 

PRINT 

6360 

PRINT 

6370 

PRINT 

6380 

PRINT 

6390 

PRINT 

64  00 

PRINT 

6410 

PRINT 

6420 

RETURN 

1 0000 

REN 

1 00 1 0 

REN  * 

1 0020 

REN 

10050 

HOME 

10060 

PRINT 

"ANN.  OS<M  COSTS=  *"; 
"FIRST  YEAR  SAVINGS^ 


YEARS" 

»"?    INT  (NO 


INT  <0M) 

*";    INT  (81) 

"DISCOUNTED  PAYBACK*    ";    INT  (DP)" 

"NET   SAVINGS    IN  CURRENT  DOLLARS* 

"ASSUMED  DISCOUNT  RATE=    ";(D  +   .08)    *    100 "7." 

"COST  OP  ELECTRICITY*  *"CG"/KWH" 

"FUEL  PRICE  ESCALATION  RATE* 

"INFLATION  RATE*  8'/." 
D$; "PR#0" 

GET  AND  READ  CITY  FILE  * 


(ES   *    100)    +  8'"/." 


INVERSE  :  GOSUB  21380 
*  ENTER  THE  NUMBER  OF 
AND  PRESS  RETURN   "; i 


1 0070 
10090 
10100 

10110 


R  LOCATION 
:    VTAB  6 
A  =   lsB  *   10: LM 
FOR   I   =  A  TO  B: 
IF  B  <   30  THEN  A  = 
HOME   :    GOTO  10090 


THE  CITY  NEAREST   ":    PRINT   "  YOU 
GOSUB  21380:    GOSUB  21380:  NORMAL 


0:TM  =  5:BM  *  20: WW  =  13:  GOSUB 
PRINT  I; 


A 


";FI*(I):  NEXT 


10:  B 


TEXT   :    VTAB  20: LO  =  1 


22070 

GOSUB  22120 
HI   =  30:    GOSUB  20170:    REM  CHOOS 


B  +    10: LM  *  LM  +  13 


101  15 
10120 
1 0 1 30 
1 0 1 50 
10160 
10170 
10180 
10190 
1 0390 
1 0400 
1 04 1 0 


GOSUB  21140: 
I 


REM     OPEN  AND  READ 


NEXT 


PRINT 
CITY 
RE*  =  FI*(A) 
GOSUB  21120: 
INPUT  I 
FOR  J  =   1  TO 
:    INPUT  A  (J) 
NEXT  J 
GOSUB  21160 
RETURN 

REM     *    INITIALIZE  CITY 
FOR   I   «   1   TO  30:    READ  FI*(I): 

DATA       ANACONDA, BILLINGS, BOZEMAN, BROWNING, BUTTE, CHOTEAU, C0LS1  h 
IFyDILLON,ENNIS,FT  BENTON 
10420     DATA     GLASGOW ,  GLEND I VE ,  GT  FALLS,  HAMILTON,  HARLOW"!  0(4, 1  irsVk:.  ,  •  '  i 

A,  JORDAN,  KAL.  I  SPELL  ,  LEWI  ST  OWN 
10430     DATA         LIBBY,  LIVINSSTN,  MILES  C- !  >  ,!  :j  Si  UiJU -!  L  1  :  ■ 
LODGE ,  S I DNEY ,  TPSN  F  Lb ,  I-  .,  i    £  1  ONL 
RETURN 
REM 

REM     *  READ  DEFAULTS  * 
REM 

READ  A,  ZP,  PY,  TH,  PA,  A  A ,  OE ,  Z  ,  BQ ,  S  A ,  X  ( 1  )  ,  Y  ( 2 )  ,  CG,  TX 
17,4,2,60,2,2,2, 1,2,7.43225039, .68, .72, 
"OWNER  BUILT  AND  INSTALLED" 
"COMMERCIAL  PURCHASE  AND  INSTALLATION 
"COMMERCIAL  PURCHASE, OWNER  INSTALLATION" 
" BANK  LOAN  OR  OTHER  COMMERCIAL  LOAN" 
"ALTERNATIVE  INVESTMENT" 
"LOW   INTEREST  UTILITY  LOAN" 
=  "INSULATED" 
=  "UNINSULATED" 


10440 
1  0500 
10510 
10520 
1 0530 
10540 
1 0550 
1 0555 
1 0560 
10565 
10570 
10575 
10580 
1 0585 


DATA 
FA$ ( 1 ) 
FA* (2) 
FA$  <  3 ) 
AA$ ( 1 ) 
AA$ (2) 
AA$ (3) 
Z*  <  1  > 
Z*  (2) 


036, 40 


11000  VTAB  16:  PRINT  :  GOSUB  24144:  GOSUB  24202:  PRINT  :  REM  RTN= 
•  ACCEF'T 

11010     GOSUB  24142:    GOSUB  24224:  PRINT   :   REM  ESC=RESTART 

11020     GOSUB  24118:   GOSUB  24228:  PRINT  :   REM  CTRL-N=SAME  VALUES,  NEW 
CITY 

11030     GOSUB  24126:    GOSUB  24226:  PRINT   :   REM     CTRL-R=RETURN  TO  BASE  C 

ASE 

11040     GOSUB  24124:    GOSUB  24222:  PRINT   :    RETURN   :    REM  CTRL-0=0U1T 


201 00 
20 1  1  0 
20 1 20 
20 1 30 
20140 
2'.)]  50 
20160 
20170 
20 1 80 
20 1  90 
202]  0 
20220 
20230 
20240 
20250 
20370 
20380 
20390 
20400 
20410 
20420 
20430 
20440 
20450 
20460 
20470 

3 

20480 
20490 
20500 
20510 
20520 
21 100 
21  120 
21 130 
21 140 
2 1  1 50 
21 160 

21300 
21310 
21320 
21330 
21340 
21350 
21360 
21370 


X  *  >!>>  X  X  **********  * 
•(SELECT  &  GET  ** 
********  XX  X  *  *  X  *  *  * 


!i£l«1  4S  UTILITY  SUBROUTINES  * 
REM  tl  44 «  44  44  44  #44  44  44  41  44  44  ti  4-  44  fl  44  44  #### 
REM 
RE  M 
REM 
Rt.M 
REM 
REM 
REM 

PRINT  :  PRINT 
GOSUB  20370: A 
RETURN 

PRINT  :    PRINT  "ANY 
PRINT   "ANY  MORE?  "; 
PRINT   "CHANGE  WHICH 
LO  OR  A   >  HI 

RETURN 
A  =  0:L  -  Ol 
GOTO  21470 


YOUR  SELECTION  =>  "; 
VAL    (At) :    IF  A  <  LO  OR  A 


>  HI   THEN  20170 


IF  A  < 
PRINT 
A*  =  "" 
ONERR 
GET  AE* 
IF  AE* 


CHANGES 
:  GOSUB 
ONE?  ==>  " 
THEN  20230 


";:    GOSUB  20790:  RETURN 
0790:  RETURN 
";:   GOSUB  20370: A  =  VAL 


(A*) 


REM       L  -STRING  LENGTH 


AE* 
AE* 
AE* 
AE* 


IF 
IF 
IF 
IF 
REM 
REM 
IF  AE* 

:A*  = 
PRINT 
REM 

A*  =  A* 
IF  L  - 
RETURN 
FOR  I 
PRINT 
PRINT 
PR  I  NT 
PR  I  NT 
PRINT 

GOSUB 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
HOME 


CHR* 
CHR* 
CHR* 
CHR* 
CHR* 


(13) 
(27) 
(14) 
(18) 
(17) 


THEN 
THEN 
THEN 
THEN 
THEN 


PRINT  I 

RETURN 

RETURN 

RETURN 

END 


RETURN 


LEFT* 
AE*; 


CHR*    (8)  THEN 
(A*,L):  GOTO 


PRINT 

20390 


CHR*  (8) 


GOSUB  20510: L  =  L  - 


+  AE*:L   =  L  +   1:   GOTO  20390 

1   <     =0  THEN     POP  :   GOTO  20370 


=   1   TO  350:    NEXT   :  RETURN 
D*;"OPEN   ";FI*(A>:  RETURN 
B$ ; "DELETE   ";FI*(A>:  RETURN 
D* ; "READ   ";FI*(A):  RETURN 
D*; "WRITE   ";FI*(A):  RETURN 
D*; "CLOSE   "FI*(A>:    PRINT  : 
1100:  RETURN 


PRINT     TAB (    10) "FILE  CLOSED' 


******** ********** 

*  MODULE  HEADERS  * 
****************** 


HTAB  4:  INVERSE 


PRINT  HDR* (HDR) :  NORMAL 


PRINT 


RETURN 


21380     GOSUB  22200:   VTAB  Y5  +   1 :    HTAB  X5  +   1 :   FOR  1=1   TO  39  -  X5:  FRINT 

"   ";:   NEXT   :   PRINT   :  RETURN 
21400  REM 
21410  REM 

21420  REM  ****************** 
21430  REM  *  ERROR  HANDLING  * 
21440     REM  ****************** 

21450  REM 
21460  REM 

21470     REM  *     ERROR  HANDLING, DOS  * 

21480  X  =     PEEK    (222)  :    POKE  216,0 

21500     IF  X  =  7  THEN     PRINT   "VOLUME  MISMATCH  IN  LINE   ";   PEEK   (218)  + 

256  *     PEEK   <219);"   **":  STOP 
21502     IF   X  =   1   THEN     PRINT    "LANGUAGE  NOT  AVAILABLE":  STOP 
21504     IF  X  ■=  2  OR   X  =  3  THEN     PRINT   "RANGE  ERROR   IN  LINE   ";    PEEK  (21 

8)   +  256  *     PEEK    (219);"   **":  STOP 
21506     IF  X  =  4  THEN     PRINT    "DISK  WRITE  PROTECTED":  STOP 

21508     IF   X  =  5  THEN     PRINT    "END  OF  DATA   IN  LINE   ";    PEEK    (218)    +  256  * 

PEEK   (219);"   **":  STOP 
21510     IF  X  =  6  THEN     PRINT   "FILE  NOT  FOUND  IN  LINE   ";   PEEK   (218)   +  2 

56  *     PEEK   (219);"   **":  STOP 
21512     IF  X  =  8  THEN     PRINT    "I/O  ERROR":  STOP 
21514     IF  X  =  9  THEN     PRINT   "DISK  FULL":  STOP 
21516     IF   X  =   10  THEN     PRINT    "FILE  LOCKED":  STOP 

21518     IF  X  =   11   THEN     PRINT    "SYNTAX   ERROR   IN  LINE          PEEK    (218)    +  25 

6  *     PEEK   (219);"   **":  STOP 
21520     IF  X  =   12  THEN     PRINT   "NO  BUFFERS  AVAILABLE   IN  LINE         PEEK  (2 

18)   +  256  *     PEEK    (219);"   **":  STOP 
21522     IF  X  =   13  THEN     PRINT   "FILE  TYPE  MISMATCH  IN  LINE   ";   PEEK  (218 

)    +  256  *     PEEK    (219);"   **":  STOP 
21524     IF  X  =   14  THEN     PRINT   "PROGRAM  TOO  LARGE":  STOP 
21526     IF  X  =   15  THEN     PRINT    "NOT  DIRECT  COMMAND":  STOP 
21528     IF  X  =   16  THEN     PRINT   "SYNTAX  ERROR   IN  LINE   " ;   PEEK   (218;   +  25 

6  *     PEEK   (219);"   **":  STOP 
21530     IF   X  =  22  THEN     PRINT    "RETURN  WITHOUT  GOSUB   IN  LINE   ";    PEEK  (2 

IB)   +  256  *     PEEK    (219);"   **":  STOP 
21532     IF   X  =  42  THEN     PRINT    "OUT  OF  DATA   IN  LINE   ";    PEEK    (218)    +  256 

*     PEEK    (219);"    **":  STOP 
21534     IF   X  =  53  THEN     PRINT    "ILLEGAL  QUANTITY   IN  LINE   ";    PEEK    (218)  + 

256  *     PEEK    (219);"   **":  STOP 
21536     IF   X  -  69  THEN     PRINT   "OVERFLOW  ERROR   IN  LINE   ";    PEEK    (218)  + 

256  *     PEEK    (219);"   **":  STOP 
21538     IF  X  =  77  THEN     PRINT   "OUT  OF  MEMORY  IN  LINE   ";   PEEK   (218)   +  2 

56  *     PEEK   (219);"   **":  STOP 
21540     IF  X  =  90  THEN     PRINT   "UNDEFINED  STATEMENT   IN  LINE   ";   PEEK  (21 

8)   +  256  *     PEEK    (219);"   **":  STOP 
21542     IF  X  =   107  THEN     PRINT   "BAD  SUBSCRIPT   IN  LINE  ";   PEEK   (218)  + 

256  *     PEEK    (219);"    **":  STOP 
21544      IF   X   =   120  THEN     PRINT   "RED I MENS I ONED  ARRAY   IN  LINE   ";    PEEK  (2 

18)   +  256  *     PEEK    (219);"   **":  STOP 
21546     IF   X   =   133  THEN     PRINT   "DIVISION  BY   ZERO   IN  LINE    ";    PEEK  (218) 

+  256  *     PEEK    (219);"    **":  STOP 
21548     IF   X  =   163  THEN     PRINT    "TYPE  MISMATCH   IN  LINE   ";    PEEK    (218)  •+ 

256  *     PEEK    (219);"    **":  STOP 
21550     IF   X  =   176  THEN     PRINT    "STRING  TOO  LONG   IN  LINE   ";    PEEK    (218)  + 

256  *     PEEK    (219);"    **":  STOP 
21552     IF   X  =   191   THEN     PRINT    "FORMULA  TOO  COMPLEX    IN  LINE   ";    PEEK  (2 

18)    +  256  *     PEEK    (219);"    **":  STOP 
21554      IF   X   =  224  THEN     PRINT    "UNDEFINED  FUNCTION   IN  LINE   ";    PEEK  (21 
8)    +  256  *     PEEK    (219);"    **":  STOP 


21556 

2155B 

21560 
2 1 800 
21B10 
21820 
21830 
21840 
2 1 850 
21860 
21870 
21880 
21890 
22000 
220 1 0 
22020 
22030 
22040 
22050 
22060 
22070 
22080 
22090 
22100 
221 10 
22 1 20 
22130 
22140 
22160 
22170 
22 1 80 
22200 
24000 
24002 
24004 
24006 
24008 
24010 
24012 
24014 

24016 

24018 

24100 
24  1  02 
24  1  04 
24  1  06 
24  1  OS 
24 1 1  0 
241  12 
241  14 
241  16 
241  18 
24 1  20 
24  1 22 
24124 


***************  ****** 

*  PRINTER  AND  DRIVE  * 
********************* 

"DO  YOU  WANT  A  HARD  COPY? 
-   "N"   THEN  RETURN 
D*;"PR#"   +     STR*  (PS): 
D*;"PR#0":  RETURN 


***************** 

*  SCREEN  SETUPS  * 
***************** 


34, TM: 


34,  TM: 

35,  BM: 

32,  LM: 

33,  WW: 


RETURN  :  REM 
RETURN  :  REM 
POKE  33,40 


IF  X  *  254   THEN  PRINT 
218)    +  256   *     PEEK  (219) 

IF  X  -  255  THEN  PRINT 
PEEK    (219);"    **":  STOP 
GOTO  980 
REM 
REM 
REM 
REM 
REM 
REM 
PRINT 
IF  A* 
PRINT 
PRINT 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
POKE 
REM 
REM 
POKE 
POKE 
POKE 
POKE 
REM 
REM 
REM 
POKE 
X5  = 
REM 
REM 
REM 
REM 
REM 
REM 
REM 

INVERSE 
CHR* (13) 

INVERSE 
CHR* (32) 

INVERSE 

REM     CHR* (27) 

INVERSE 

INVERSE 

INVERSE 

INVERSE 

INVERSE 

INVERSE 

INVERSE 

INVERSE 

I N VERSE 

INVERSE 

INVERSE 

INVERSE 

INVERSE 


"BAD  RESPONSE  TO  INPUT   IN  LINE   ";  PEEK 
;"   **":  STOP 
"CONTROL-C   IN  LINE 


PEEK    (218)  + 


256 


;   GDSUB  20790 
RETURN 


POKE  35, BM:    POKE  32, LM:    POKE  33, WW:  RETURN 


TOP 

BOTTOM 
LM:  RETURN 


RETURN   :    REM  WINDOW  WIDTH 


REM 


LEFT  AND  WIDTH 


36, HC:  RETURN 
PEEK    (36) : Y5  • 


REM  HORIZONTAL  CURSOR 
PEEK 


(37):  RETURN 


********************** 
*  SINGLE  KEY  OPTIONS  * 
********************** 


PRINT 
PRINT 
PR  I  NT 


"RTN" ; 
"SPC"; 
"ESC": 


NORMAL 
NORMAL 
NORMAL 


PRINT 
PRINT 


"= ACCEPT 
"=NEXT  EVT 


:    RETURN   :  REM 
":    RETURN   :  REM 


PRINT    "=NEXT  SWMR 


RETURN 


PRINT 

"CTRL 

-A"  : 

NORMAL  : 

RETURN 

REM 

CHR* ( 1 ) 

H 

PRINT 

"CTRL 

-B" 

NORMAL 

•  RETURN 

REM 

CHR* (2) 

PRINT 

"CTRL 

-C"  : 

■ 

NORMAL  i 

RETURN 

* 
■ 

REM 

CHR* (3) 

PRINT 

"CTRL 

-D" 

NORMAL 

RETURN 

REM 

CHR* (4) 

PRINT 

"CTRL 

-E"  : 

NORMAL  : 

RETURN 

■ 

REM 

CHR* ( 5 ) 

PRINT 

"CTRL 

-F" 

■ 

NORMAL 

RETURN 

m 
■ 

REM 

CHR* (6) 

PRINT 

"CTRL 

-I  ": 

NORMAL  : 

RETURN 

REM 

CHR* ( 9 ) 

■ 

PRINT 

"CTRL 

-K": 

NORMAL  : 

RETURN 

• 

REM 

CHR* (11) 

PRINT 

"CTRL 

-L  "  : 

■ 

NORMAL  : 

RETURN 

■ 

REM 

CHR* (12) 

PRINT 

"CTRL 

-N  "  ; 

NORMAL  : 

RETURN 

REM 

CHR* (14) 

PRINT 

"CTRL- 

-0"; 

NORMAL  : 

RETURN 

■ 

REM 

CHR*(15) 

PRINT 

"CTRL 

-P" ; 

• 

NORMAL  : 

RETURN 

m 
m 

REM 

CHR* (16) 

PRINT 

"CTRL- 

-Q"i 

■ 

■ 

NORMAL  : 

RETURN 

: 

REM 

CHR* (17) 

24126 

24128 

24130 

24132 

24 1 34 

24136 

24 1 38 

24140 

24142 

24144 

24146 

24148 

24 1 50 

24152 

24200 

24202 

24204 

24206 

24208 

24210 

24212 

24214 

24216 

24218 

24220 

24222 

24224 

24226 

24228 

38000 

380 1 0 

38020 

3B030 

38040 

38050 

38060 

38100 

38110 

38200 

38210 

38220 

38230 

38240 

38250 

38260 

38300 

38310 

38320 

39000 

39010 

39020 

39030 

39040 

39050 

39060 

39100 

39105 

391 10 


I NVERSE 
INVERSE 
INVERSE 
INVERSE 
INVERSE 
INVERSE 
INVERSE 
REM 

INVERSE 
INVERSE 
INVERSE 
REM 
REM 
REM 
REM 
PR  I  NT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
REM 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
BR  I  NT 
PRINT 

PRINT 
PRINT 
PRINT 


"  CTRL 

-R  " ; 

NORMAL  : 

RET!  IPN  ■ 

r\t  II 

PUD*  /  4  rt  \ 

LHK* (IB) 

"CTRL 

-S"  ; 

1 

NORMAL  ■ 

FFTI  IfvKI  - 
r\ d  i  ur\iv  : 

rvtrl 

CHR$  < 1 9) 

"CTRL 

-T" ; 

NORMAL  : 

RFTI  IRI\J  • 

nc  i  Unix  ■ 

C'CM 

ntn 

LHR*  (2U) 

"CTRL- 

-V" ; 

NORMAL  - 

htn 

LHR*  <  22 ) 

"CTRL 

-W"; 

■ 

i 

NORMAL  : 

RETURN  : 

REM 

CHR*  <23) 

"CTRL- 

-Y"; 

NORMAL  : 

RETURN  : 

REM 

CHR$(25) 

"CTRL 

-Z"; 

■ 

NORMAL  : 

RETURN  : 

REM 

CHR$(26) 

"ESC" ; : 
"RTN" ; : 
"SPC"; : 


NORMAL 
NORMAL 
NORMAL 


RETURN 
RETURN 
RETURN 


REM 
REM 
REM 


CHR$<32> 
CHR*(13) 
CHR*<32) 


"=ACCEPT  VALUE   ";:  RETURN 
"=DEFAULT  VALUE   ";:  RETURN 
"=NEXT   ITEM   ";:  RETURN 
"= INSERT    ";:  RETURN 
"=DELETE    ";:  RETURN 
"=ADD   ";:  RETURN 
"=REPLACE   ";:  RETURN 
"=CONTINUE   "; :  RETURN 
"=CALCULATE    ";:  RETURN 
"=RETURN  TO  MENU   ";:  RETURN 
"=QUIT";:  RETURN 

"=RESTART, VALUES   INTACT   ";:  RETURN 
"^RETURN  VALUES  TO  BASE  CASE    ";:  RETURN 
"=SAME  VALUES,    NEW  CITY";:  RETURN 


*************** 

*  GET  VALUE  1  * 
*************** 


22200:  REM  XI, Yl 
CV:  RETURN 


*************** 
*  GET  VALUE  2  * 
*************** 


REM 
REM 
REM 
REM 
REM 
REM 
GOSUB 
PRINT 
REM 
REM 
REM 
REM 
REM 
REM 
REM 

VTAB  Y5  +  1: 
GOSUB  20370: 
RETURN 
REM 
REM 
REM 
REM 
REM 
REM 
REM 

IF  A*  < 
VTAB  Y5  + 
IF  AE$  = 
AND  A*  < 


HTAB  X5  +  1 
GOSUB  39100 


********************* 
*  KEY  TESTS,  EVENTS  * 
********************* 


""   THEN  LD  = 
1:   HTAB  X5  +  1 
CHR$    (13)  AND 
>   ""   THEN  CV  ■ 


39115     IF  A*  < 


LEN    (  STR* 
+  LD:  CALL 
VAL    <  A* )  > 
VAL  (A*) 


<    VAL  (A*))) 

-  868 

-  LO  AND  VAL 


(A*)  < 


H 


AND     VAL    <A$)    <   LO  OR     VAL    (A*)    >  HI   THEN  POP 


GOTO  38300 

39120      IF  AE*  =     CHR*    (27)    THEN     HOME   :    POP   :    POP   :    GOTO  680 
39130      IF  AE*  =     CHR*    (14)    THEN     GOSUB   10050:    POP   :    POP   :    GOTO  680 
39140      IF  AE*  =     CHR*    (18)    THEN     RESTORE    :    GOSUB   10400:    GOSUB   10530s  POP 

:    POP   :    GOTO  680 
39150  RETURN 

39300     PRINT   "        CURRENT  VALUE  m>   " ; :    GOSUB  22200:    PRINT  CV:LD  -  LEN 
(   STR*    (CV)>:    GOSUB   11000:    REM       KEY  TESTS 


"^9"^  1  O 

GOSUB 

38300: 

RETURN 

4  pppp 

RETURN 

r~  ;")/")  i')/') 
\J\J\f\J\J 

REM 

uJ      v*  1  v 

REM 

5  O  O  2  (j 

REM  J***************** 

u 

REM  * 

GOSUB 

LISTING  * 

50040 

REM  ***#****#•**♦**♦*♦ 

50050 

RFM 

PPM 

wJ  V  K-'  /  \.' 

T>M  1  7f')i 

PPM 

YUUK   b^Ltu 1  ION? 

t j  \_i  O  LJ  £? 

^_  U  ^_  1 s 

PPM 

ANY  CHANGES? 

5  Of)  9  m 

GOSUB 

r\ui  1 

MIN  I     1  lurvt  : 

50 1 00 

GOSUB 

20230  ■ 

R'FM 

"PWAMRP     lilUITf^UJ     nhlC^  11 
l^rtHINOt    WnlLM    UrJt  : 

501 1 0 

REM 

501  50 

REM 

SO  1  AO 

JL O  /  - 

PPM 

rtfc.Pi 

Ml  II  TT  C'l  c    r~  l — r 

rlUL  1  J.  1-  I_L  GET 

50 1 70 

iUo  /0 : 

REM 

ANY  KEY  TO  CONTINUE 

SO  1  OO 

RDCil  IR 

L3U3LJC* 

Tl7QA  • 

C'PM 

r\h.rl 

II  V/  /  Kin  11 

Y  /  N  r 

50 1 85 

^Ud/iU  : 

I S  TH I S  OK? 

50 1 90 

ppm 

50200 

RDQI  IR 

LJlXOUC 

51 1  riAf">  ■ 
1  : 

PPM 

PA      C"  T  1    C      nCifM  1 

KM   I-  I  Lb  OPEN 

5021 0 

ROQI  IR 

*?  1  H7ri . 
-t-  ±  k>  /  : 

C'PM 

P/\      PTI  1 —  r— ,  1 —  y-\  y-T, 

KA   F-ILE  READ 

'  a. 

nsncii  ir 

L3  L_»  O  LI 

^  x  uou : 

C'PM 

KA  FILE  WRITE 

i_i ' i 

oUbUB 

21090: 

REM 

RA  FILE  CLOSE 

50240 

r\c.l  1 

.  1 '  1  .•  .  n  t 

CflQI  IR 

1 1  w : 

C'CM 

T~\  r —  1    /\\/      r—  r~i  r~i      i-v  t  r-v  r-.  1     a  \  ; 

DELAY  FOR  DISPLAY 

PPM 
rv  t  II 

50270 

rnni  it-, 

211 20 : 

REM 

SEQUENTIAL  FILE  OPEN 

50280 

p,nci  id 

^11 jU : 

C'CM 

REM 

SEQUENTIAL   FILE  DELETE 

pnci  itj 
bubUb 

1 40 : 

REM 

SEQUENTIAL  FILE  READ 

/"I  T  i~\  t~\ 

nnQi  ir 

REM 

SEQUENTIAL  FILE  WRITE 

bUOUtJ 

ill oO! 

C'CrM 

Ktn 

SEQUENTIAL   FILE  CLOSE 

50320 

n.  t(  1 

50/^/^  O 

RflQI  IR 
CJlJb>Uci 

^.1 1  /  0 : 

C'CM 

M0N   C , I , 0 

«J  U  ■_■  *T  <  J 

bUbUo 

ill  00 : 

r  ■  j —  ^ 

REM 

N0M0N  C , I , 0 

50350 

RFM 

SO"?  A/") 
*-J 0 

GOSUB 

21360: 

REM 

PRINT  MODULE  HEADER 

SO  "^71") 

GOSUB 

21460: 

REM 

0NERR  GOTO,  DOS 

SOXQO 

GOSUB 

21730: 

REM 

CONVERT  TIME  TO  DECIMAL 

S  O  A.  (~)  ("> 

REM 

504  1  O 

GOSUB 

21850: 

REM 

"HARD  COPY?" 

GOSUB 

21880: 

REM 

TURN  ON  PRINTER 

50430 

GOSUB 

21890: 

REM 

TURN  OFF  PRINTER 

50440 

REM 

50450 

GOSUB 

22070: 

REM 

ALL  SCREEN  POKES 

50460 

GOSUB 

22 1 00 : 

REM 

SET  TOP  MARGIN 

50470 

GOSUB 

221  ]0: 

REM 

SET  BOTTOM  MARGIN 

50480 

GOSUB 

22120: 

REM 

SET  LEFT  MARGIN 

50490 

GOSUB 

22130: 

REM 

SET  WINDOW  WIDTH 

50500 

GOSUB 

22150: 

REM 

SET  CURSOR,    H0R  AND  VERT 

50510 

GOSUB 

22180: 

REM 

SET  CURSOR,  H0R 

50520 

GOSUB 

22190: 

REM 

SET  CURSOR,  VERT 

50530 

GDSUB 

24014: 

50540 

GOSUB 

24016: 

50550 

GOSUB 

24018: 

50560 

GOSUB 

24 1 00 : 

50570 

GOSUB 

24102: 

50580 

GOSUB 

24 1 04 : 

50590 

GOSUB 

24 1 06 : 

50600 

GOSUB 

24108: 

50610 

GOSUB 

241 10: 

50620 

GOSUB 

241 12: 

50630 

GOSUB 

241 14: 

50640 

GOSUB 

24116: 

50650 

GOSUB 

241 18: 

50660 

GOSUB 

24120: 

50670 

GOSUB 

24122: 

50680 

GOSUB 

24124: 

50690 

GOSUB 

24126: 

50700 

GOSUB 

24123: 

50710 

GOSUB 

24130: 

50720 

GOSUB 

24132: 

50730 

GOSUB 

24134: 

50740 

GOSUB 

24136: 

50750 

GOSUB 

24138: 

50760 

REM 

50770 

GOSUB 

24 142: 

50780 

GOSUB 

24144: 

50790 

GOSUB 

24146: 

50800 

REM 

50810 

GOSUB 

242Q2 • 

50820 

GOSUB 

24204: 

50830 

GOSUB 

24206: 

50840 

GOSUB 

24208: 

50850 

GOSUB 

24210: 

50860 

GOSUB 

24212: 

50870 

GOSUB 

24214: 

50880 

GOSUB 

24216: 

50S90 

GOSUB 

24218: 

50900 

GOSUB 

24220: 

50910 

GOSUB 

2  A  22" 

5 1 000 

GOSUB 

24014: 

51010 

REM 

REM 

"RTN= ACCEPT " 

REM 

"SPONEXT  EVT" 

REM 

"Esc=ne; 

<T  SWMR 

REM 

CTRL-A 

1 

REM 

CTRL-B 

2 

REM 

CTRL-C 

3 

REM 

CTRL-D 

4 

REM 

CTRL-E 

5 

REM 

CTRL-F 

6 

REM 

CTRL- I 

9 

REM 

CTRL— K 

1 1 

REM 

CTRL-L 

12 

REM 

CTRL-N 

14 

REM 

CTRL-0 

15 

REM 

CTRL-P 

16 

REM 

CTRL-Q 

17 

REM 

CTRL-R 

18 

REM 

CTRL-S 

19 

REM 

CTRL-T 

REM 

CTRL-V 

REM 

CTRL-W 

X-  ■_> 

REM 

CTRL-Y 

25 

REM 

CTRL-Z 

26 

REM 

ESC  32 

REM 

RTN  13 

REM 

SPC  32 

REM  "ACCEPT" 

REM  "DEFAULT" 

REM  "NEXT  ITEM" 

REM  " INSERT" 

REM  "DELETE" 

REM  "ADD" 

REM  "REPLACE" 

REM  "CONTINUE" 

REM  "CALCULATE" 

REM  "RTN  TO  MENU" 

REM  "QUIT" 

REM  "RTN= ACCEPT; ESC=QUIT 


15  copies  of  this  public  document  were  published 
at  an  estimated  cost  of  $5.00  per  copy,  for  a  total 
cost  of  $75.00,  which  includes  $75.00  for  printing 
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